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ABSTRACT 


.The  general  problem  of  the  interaction  between  a  mon¬ 
atomic  gas  and  a  solid  surface  is  investigated  from  a  mathe¬ 
matical  point  of  view  by  the  use  of  a  scattering  kernel. 

This  general  scattering  kernel  must  fulfill  the  obvious  con¬ 
ditions  of  normalization  and  non-negativity  and  has  recently 
been  shown  to  also  satisfy  the  condition  of  detailed  balanc¬ 
ing  or  reciprocity.  Previous  gas-surface  theories,  both  ana¬ 
lytical  and  empirical,  are  considered  in  view  of  reciprocity. 
It  is  shown  that  most  of  these  descriptions  are  satisfactory 
for  gross  gas-surface  interactions  but  fail  to  describe  the 
detailed  interaction.  By  assuming  a  .product  solution  in  each 
of  three  coordinate  systems,  an  infinite  series  of  solutions 
is  obtained.  The  first  scattering  kernel  in  each  series  is 
obtained  in  closed  mathematical  form,  while  the  second  and 
third  are  evaluated  numerically.  Recurrence  relations  are 
obtained  for  the  coefficients  in  the  infinite  series  for  all 
solutions.  It  is  shown  that  the  first  scattering  kernel  in 
rectangular  coordinates  adequately  describes  the  experiment¬ 
ally  observed  results  of  gas-surface  interaction.  This  solu¬ 
tion  has  two  parameters  which  describe  the  gas^surface  inter¬ 
action.  It  is  shown  that  these  parameters  are  related  to  a 
tangential  and  a  normal  thermal  accommodation  coefficient. 

The  scattering  kernels  obtained  are  integrated,  in  some  cases 
analytically  and  others  numerically,  with  various  weighting 
functions  in  order  to  obtain  mean  reflected  properties,  both 
global  and  local,  such  as  molecular  flux,  speed,  and  energy. 
These  various  results  are  presented  for  typical  gas-surface 
interaction  parameters.  It  is  shown  that  some  combination  of 
parameters  describe  general  scattering  characteristics.  The 
theoretical  results  are  in  agreement  with  existing  detailed 
measurements  of  reflected  quantities  such  as  intensities  and 
mean  speed,  for  a  wide  range  of  incident  energy  and  direction. 
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I.  INTRODUCTION 

The  interaction  of  a  gas  with  a  solid  surface  is  very 
important  in  several  fields  of  science.  The  aerodynamicist 
would  like  to  know  the  gross  interaction  effects  such  as  the 
total  momentum  and  energy  exchange  under  conditions  such  as 
the  flight  of  an  artificial  satellite  through  the  rarefied 
upper  atmosphere.  The  mathematician,  working  in  the  field 
of  kinetic  theory,  desires  a  mathematical  description  of  the 
gas-surface  interaction  for  a  boundary  condition  to  the 
Boltzmann  equation.  Theoretical  and  experimental  results  on 
gas-surface  interaction  help  the  physical  chemist  better 
understand  the  complex  process  of  surface  catalyticity  and 
help  the  physicist  understand  the  surface  structure  and  the 
gas-solid  interaction  potentials. 

In  general,  the  knowledge  and  understanding  of  gas- 
surface  interaction  phenomena  are  less  advanced  than  that  in 
the  bulk  gas,  this  in  spite  of  the  fact  that  most  experiment¬ 
ation  in  the  bulk  gas  is  influenced  to  a  large  extent  by 
solid  boundaries.  It  is,  in  general,  true  that  the  knowledge 
of  the  interface  phenomena  between  two  disciplines  proceeds 
slower  than  that  in  either.  Usually,  it  is  necessary  to 
resort  to  empirical  methods  for  descriptions  of  these  phenom¬ 
ena,  e.g.,  heat  transfer  coefficients. 

Many  times  in  this  paper  and  other  papers  the  terms 
"reflected  gas"  and  "reemitted  gas"  are  used  interchangeably 
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to  mean  the  gas  after  it  has  interacted  with  the  surface, 
irrespective  of  the  type  of  interaction. 

1.1  HISTORICAL  BACKGROUND 

It  is  now  approaching  100  years  since  Knudt  and  War¬ 
burg  [l]1  in  1875  proposed,  based  upon  measurements  of  the 
damping  of  a  vibrating  disk  in  a  rarefied  gas,  that  under 
certain  conditions  a  gas  may  slip  adjacent  to  a  solid  bound¬ 
ary.  They  speculated  that  something  analogous  may  occur  in 
rarefied  gases  when  heat  is  being  conducted,  namely,  a 
temperature  jump  at  a  solid  boundary.  Prior  to  this  time  it 
had  been  assumed  that  a  gas  in  contact  with  a  solid  surface 
has  identical  velocity  and  temperature  as  the  surface,  i.e., 
the  zero  slip  condition.  Osborne  Reynolds  [2]  in  1879  intro¬ 
duced  a  factor  of  slip  between  mean  tangential  velocities  of 
incident  and  reflected  molecules.  He  recognized  that  this 
factor  somehow  depended  upon  the  nature  of  the  gas-surface 
interaction.  He  also  realized  that  incomplete  molecular  ac¬ 
commodation  is  not  a  prerequisite  for  velocity  slip  or  tem¬ 
perature  jump  at  a  surface.  However,  this  factor,  the  in¬ 
verse  of  which  he  considered  a  coefficient  of  roughness,  con¬ 
cerns  macroscopic  properties. 

Maxwell  [3]  later  that  same  year  was  the  first  to 
propose  a  gas  surface  interaction  model  based  upon  the  kinetic 
theory  of  gases,  which  he  pioneered  [4].  He  showed  that  both 

lumbers  in  brackets  refer  to  similarly  numbered 
references  in  the  List  of  References. 
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a  perfectly  elastic  smooth  surface  and  a  surface  of  such  pro¬ 
perties  that  the  gas  being  reflected  has  the  same  distribu¬ 
tion  of  velocities  as  if  it  were  emitted  from  a  static  gas  at 
the  surface  temperature  (complete  accommodation)  do  not  dis¬ 
turb  a  gas  in  equilibrium  with  the  surface.  He  proposed: 

...  to  treat  the  surface  as  something  intermedi¬ 
ate  between  a  perfectly  reflecting  and  a  perfectly 
absorbing  surface,  and,  in  particular,  to  suppose 
that  of  every  unit  of  area  a  portion  f  absorbs  all 
the  incident  molecules,  and  afterwards  allows  them 
to  evaporate  with  velocities  corresponding  to  those 
in  still  gas  at  the  temperature  of  the  solid,  while 
a  portion  1-f  perfectly  reflects  all  the  molecules 
incident  upon  it. 

Even  today  the  gas-surface  model  of  Maxwell  is  pro¬ 
bably  the  most  widely  used  for  calculating  gross  gas-surface 
effects.  It  is  interesting  that  his  original  paper  [3]  did 
not  contain  any  mention  of  the  gas-surface  interaction;  this 
was  added  in  an  appendix  at  the  insistence  of  a  referee. 

Smoluchowski  [5]  was  the  first  to  measure  the  heat 
transfer  between  rarefied  gases  and  a  surface  and  to  show 
that  there  is  indeed  a  temperature  jump.  He  also  recog¬ 
nized  the  need  for  parameters  to  indicate  the  degree  that 
gas  molecules,  upon  an  encounter  with  a  surface,  assume 
the  conditions  of  the  surface.  These  parameters,  which 
can  be  defined  for  each  molecular  quantity  that  can  be  af¬ 
fected  by  the  encounter,  are  known  as  accommodation  coef¬ 
ficients.  However,  the  thermal  accommodation  coefficient 
as  defined  by  him,  was  awkward  since  it  could  vary  over  all 
positive  numbers.  His  unlucky  choice  of  the  definition  of 
this  factor,  caused  Knudsen  [6]  to  be  widely  recognized 
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as  the  originator  of  the  thermal  accommodation  coefficient, 
a,  which  he  defines  to  vary  between  zero  and  unity. 

Knudsen  experimentally  determined  a  for  several  different 
gas-surface  combinations  and  found  that  it  varied  from  0.26 
for  hydrogen  in  contact  with  polished  platinum  up  to  ap¬ 
proximately  0.87  for  the  heavier  gases  such  as  nitrogen. 

Both  Smoluchowski  and  Knudsen  defined  accommodation 
in  terms  of  temperature  since  they  were  especially  interested 
in  flows  close  to  equilibrium  for  which  temperature  could  be 
defined.  Any  definition  of  thermal  accommodation  involving 
temperature  implicitly  assumes  that  the  incident  and  re¬ 
flected  gas  has  a  temperature  associated  with  it,  i.e.,  it 
has  a  Maxwellian  velocity  distribution.  Recently  Kuscer, 
et  al  [7]  read  into  Knudsen 's  definition  of  a  the  assumption 
that  for  a  gas  with  a  Maxwellian  distribution  incident  upon 
a  surface,  that  the  reflected  distribution  is  also  Maxwellian 
at  a  temperature  determined  from  the  definition  of  the 
accommodation  coefficient.  This  hypothesis,  which  Kuscer 
[7]  refers  to  as  "Knudsen ’s  assumption"  will  be  compared  to 
theoretical  results  later.  However,  Knudsen  himself  recog¬ 
nized  that  for  departures  from  equilibrium  the  reflected 
molecules  may  not  have  a  Maxwellian  velocity  distribution 
and  that  the  definition  of  a  must  be  written  in  terms  of 
energy. 

The  experimental  results  of  Smoluchowski  [5]  and 
Knudsen  [6]  and  also  those  of  Soddy  and  Berry  [8]  and 
Blodgett  and  Langmuir  [9]  for  the  thermal  accommodation  of 
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various  gas-surface  interactions  indicated  two  general 
trends.  First,  as  was  probably  expected,  a  was  lower  for 
polished  surfaces  than  for  rough  surfaces.  Second,  a  was 
lower  for  lighter  gases  and  approached  unity  for  heavier 
gases. 

Baule  [10]  in  1914  was  the  first  to  attempt  a  micro¬ 
scopic  analysis  of  gas-surface  interactions.  He  assumed 
that  both  the  gas  molecules  and  the  surface  atom  behave  as 
hard  spheres  with  a  mass  ratio  p,.  The  result  that  he  ob¬ 
tained  for  the  case  of  a  zero  temperature  solid  was: 

a  =  2u/(l+n)2 

This  result  was  in  qualitative  agreement  with  the  variation 
of  a  with  gas  molecular  weight  for  ^  <  1.  In  general,  how¬ 
ever,  this  simple  theory  could  not  be  expected  to  contribute 
to  a  better  understanding  of  gas-surface  interactions  since 
neither  the  surface  structure  nor  the  interaction  potential 
is  modeled  in  a  realistic  way. 

Following  Baule,  several  other  microscopic  gas- 
surface  theories  were  proposed  based  upon  both  classical 
and  the  newly  developed  quantum  mechanics.  Most  of  the 
early  analyses  were  of  necessity  one-dimensional  treatments 
which  were  either  simple  and  not  very  meaningful  or  complex 
and  not  amenable  to  mathematical  solutions.  Only  with  the 
introduction  of  large  analogue  and  digital  computers  has  it 
been  possible  to  propose  and  solve  more  realistic  three- 
dimensional  models. 
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As  the  theoretical  models  advanced  so  did  the 
experimental  capabilities.  The  experimental  results  pre¬ 
sented  in  [5,  6,  8  and  9]  were  based  upon  near  equilibrium 
conditions,  where  the  gas  incident  upon  the  surface  had  a 
near  Maxwellian  velocity  distribution  with  a  temperature 
close  to  the  surface  temperature.  Beginning  around  1930 
there  was  renewed  interest  in  experimental  gas  surface 
interactions,  which  was  prompted  by  the  prediction  of 
diffraction  effects  by  quantum  theory.  Unidirectional 
molecular  beams  were  generated  by  the  use  of  an  effusive 
orifice  or  Knudsen  source.  Most  notable  among  the  experi¬ 
mental  investigators  of  this  period  was  Estermann  [11,  12]. 

It  was  not  until  1951  that  Kantrowitz  and  Gray  [13] 
suggested  an  improvement  on  molecular  beam  generation. 

Their  aerodynamic  nozzle  and  skimmer  permitted  a  several 
fold  increase  in  molecular  flux  as  well  as  a  more  nearly 
monoenergetic  beam. 

Molecular  detection  techniques  have  also  improved  to 
the  point  such  that  today  the  velocity  distribution  function 
of  molecules  reflected  from  a  surface  can  be  determined, 
thus  permitting  the  determination  of  desired  gas-surface 
pa rameters. 

The  earlier  theoretical  gas-surface  models  were 
required  to  reproduce  only  the  gross  accommodation  coeffi¬ 
cients  of  which  a  was  the  most  important.  However,  with 
improved  experimental  techniques,  it  was  necessary  that  the 
theoretical  results  be  able  to  describe  these  more  detailed 
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experimental  data.  The  most  important  of  these  was  the 
distribution  of  reflected  molecular  flux  as  a  function  of 
reflected  direction,  usually  referred  to  as  the  "spatial 
distribution. " 

The  most  complete  classical  model  at  present  is 
probably  that  of  Oman,  et  al  [14].  This  model  requires  a 
statistical  or  Monte  Carlo  approach  in  order  to  obtain  a 
solution  and  is  therefore  not  convenient  for  a  mathematical 
description  of  the  gas-surface  interaction.  There  are  many 
other  analyses  such  as  the  "hard  cube"  model  of  Logan  and 
Stickney  [15]  and  the  "soft  cube"  model  of  Logan  and  Keck 

[16] ,  These  models  are  either  too  complicated  for  a  general 
description  of  the  gas-surface  interaction  or  do  not  ade¬ 
quately  describe  it.  For  a  more  complete  discussion  of  gas- 
surface  interaction  theories  see  the  review  paper  by  Goodman 

[17] .  Recently,  Busby  [18]  has  proposed  a  gas-surface  model 
that  utilizes  a  potential  well  at  the  surface  in  order  to 
explain  the  observed  capture  coefficient  of  a  molecular 
beam  incident  upon  a  highly  cooled  surface. 

In  order  to  represent  the  gas-surface  interaction 
in  mathematical  terms,  an  alternate  approach  has  been  to 
propose  empirical  descriptions.  As  has  been  pointed  out, 
Maxwell  [3]  proposed  a  model  of  this  sort.  More  recently, 
Nocilla  [19,  20]  has  extended  the  Maxwell  model  to  assume 
that  the  reeraitted  velocity  distribution  is  equivalent  to 
the  free  molecule  effusive  flow  of  a  gas  with  mean  velocity 
and  temperatures  to  be  determined  by  fitting  the  experimental 
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data.  Hurlbut  and  Sherman  [21]  have  applied  this  model  to 
calculate  aerodynamic  forces  on  simple  bodies  in  free 
molecule  flow.  Epstein  [22]  proposes  an  empirical  reflec¬ 
tion  model  that  matches  the  variation  of  a  with  temperature 
by  determining  three  arbitrary  constants.  These  empirical 
models  will  be  discussed  in  more  detail  later. 

It  has  been  widely  recognized  since  Maxwell  that  any 
gas-surface  model,  either  analytical  or  empirical,  must 
preserve  a  Maxwellian  velocity  distribution  at  the  surface 
temperature,  i.e.,  at  equilibrium  the  surface  does  not 
disturb  the  gas.  This  condition  has  recently  been  strength¬ 
ened  by  showing  that  not  only  is  the  Maxwellian  distribution 
maintained  at  a  surface  globally  but  also  the  principal  of 
detailed  balancing  is  satisfied.  These  conditions  will  be 
discussed  in  more  detail  in  the  next  section. 

1.2  RECIPROCITY  AND  GAS-SURFACE  INTERACTION 

The  time  independent  interaction  between  a  monatomic 
gas  and  a  solid  boundary  is  expressed  mathematically  by 
means  of  a  function  P(v’,v)  which  relates  the  probability  of 
reflected  molecular  velocity  v,  to  that  of  the  incident 
molecular  velocity,  vr,  and  the  boundary  or  wall  properties. 
Thus,  for  given  wall  conditions,  a  molecule  with  velocity  v' 

impinging  upon  the  surface  will  have  a  probability  of 

—  —  3 

P(v',v)d  v  that  it  will  be  scattered  into  the  velocity 

3  _ 

increment  d  v  about  v.  P(v%v)  is  known  as  a  scattering 
kernel  because  of  the  role  it  plays  in  Equation  3. 
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The  scattering  kernel  P(\T' ,v),  which  is  defined  only  for 
v1  directed  toward  the  surface  and  v  directed  away  is  non¬ 
negative  and  for  the  case  of  no  mass  exchange  with  the 
surface  is  normalized  in  half  velocity  space.  These  require¬ 
ments  are  expressed  mathematically  as: 

P (v' ,  v)  >  0,  *v'  *T  <  0,  v-T  >  0  (1) 

/  P(7’ ,v)  d3v  =  1  (2) 

vj>0 

See  the  List  of  Symbols  and  Figure  1,  Appendix  A,  for  the 
definition  of  symbols.  It  should  be  pointed  out  that  the 
orientation  of  the  x, z  axes  for  an  isotropic  surface  is 
arbitrary  and  is  usually  chosen  such  that  v  '  is  zero.  If 
the  surface  is  not  isotropic  but  has  differences  in  struc¬ 
ture  in  various  directions  then  the  x  and  z  axes  should  be 
defined  to  align  with  the  structure  of  the  surface  in  order 
to  simplify  the  equations. 

The  use  of  the  scattering  kernel,  as  has  been  pointed 
out  by  Shidlovskiy  [23]  is  the  most  general  representation 
of  the  boundary  conditions  at  a  solid  surface.  The  use  of 
the  scattering  kernel  is  desirable  in  order  to  completely 
specify  the  boundary  conditions  for  the  Boltzmann  equation 
when  solid  boundaries  are  involved. 
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The  velocity  distribution  function,2  f,  in  a  gas 
adjacent  to  a  solid  boundary  may  be  considered  to  consist 
of  two  parts,  those  molecules  with  velocities  directed  toward 
the  solid  surface  or  wall  and  those  with  velocities  directed 
away  from  the  wall.  Thus: 

f  =  f+  +  f" 

where:  f+  =  0,  v**j”  >  0 

f“  5  0,  <  0 

Using  the  scattering  kernel  it  can  be  shown  (Cercig- 
nani  [24],  page  51)  that  the  distribution  of  reflected  mole¬ 
cules  is  related  to  the  distribution  of  incident  molecules  by: 

f+(v)  =  -(v-jT1  /  (T’-j}p(7',"$)f“(v')  d3v'  (3) 

v-T<0 

From  consideration  of  a  gas  in  equilibrium  with  a 
surface,  it  is  obvious  that  if  f-(v*’>  is  a  static  Maxwellian 
distribution  at  the  wall  temperature,  Tw,  then  f+(v)  is  also 
a  Maxwellian  at  Tw«  This  fact  restricts  somewhat  the  func¬ 
tions  that  are  suitable  for  describing  gas-surface  interactions. 


2The  velocity  distribution  function  as  used  in  this 
paper  is  defined  as  a  probability  density  in  velocity  space 
at  a  point  in  physical  space.  This  is  the  normal  kinetic 
theory  definition  of  f.  However,  most  authors  in  the  field 
of  molecular  beams  and  gas-surface  interaction  define  a 
velocity  distribution  function  in  terms  of  flux  of  molecules 
instead  of  molecular  density.  This  frequently  causes  con¬ 
fusion  since  many  times  the  definition  used  is  not  clearly 
stated. 
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However,  from  an  equilibrium  consideration  it  can  be  shown 
that  not  only  is  the  Maxwellian  distribution  maintained  by 
the  molecular  interactions  as  a  whole,  but  each  encounter 
is  compensated  by  an  inverse  encounter.  In  other  words,  the 
general  principle  of  "detailed  balancing"  is  applicable  to 
gas-surface  interaction  as  well  as  in  the  gas  itself. 
Cercignani  [24]  recently  pointed  out  the  importance  of  de¬ 
tailed  balancing  for  the  description  of  the  gas-surface 
interaction.  Kuscer  and  Summerfield  [25]  obtained  the  same 
results  based  upon  an  analogy  with  slow  neutrons  and  referred 
to  it  as  the  "reciprocity  relation." 

This  reciprocity  relation  in  terms  of  present 
nomenclature  is: 

-  (7’  *1)  exp  (-7  2/2RT  )  P(7  ,7)  =  (7  j)  exp  (- v2/2RT  )  P(- 7,  -7  ) 

w  w 

(4) 

The  interaction  of  a  molecule  with  the  surface, 
which  is  described  by  P(7*  ,7),  takes  place  on  an  individual 
statistical  basis  and  is  independent  of  whether  the  gas  is 
in  equilibrium  with  the  wall  or  not  as  long  as  the  surface 
is  in  equilibrium  with  itself.  Therefore,  Equation  4  applies 
independent  of  the  thermodynamic  state  of  the  gas  or  the 
details  of  what  happens  on  a  microscopic  scale  during  the 
gas-surface  interaction. 

It  may  be  surprising  that  detailed  balancing  holds 
between  gas  molecules  only  when  the  gas  is  at  equilibrium 
but  is  still  applicable  for  molecular  encounters  with  the 
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surface  under  any  conditions  of  the  gas.  It  might  help 
to  think  of  detailed  balancing  as  an  interaction  between 
a  test  particle  and  an  equilibrium  reservoir.  The  test 
particle  is  a  molecule  with  velocity  v’*  before  the  en¬ 
counter  and  velocity  v  after  the  encounter.  The  reservoir 
in  one  case  is  the  equilibrium  gas  and  in  the  other  is  the 
surface.  For  non-equilibrium  conditions  in  the  gas  the 
reservoir  consisting  of  gas  molecules  ceases  to  exist  while 
the  surface  is  assumed  to  remain  unchanged.3  It  should  be 
pointed  out  that  under  varying  molecular  fluxes  to  a 
surface  that  the  characteristics  of  the  surface  may  change. 
For  example,  under  the  condition  of  high  flux,  a  relatively 
large  amount  of  gas  may  be  adsorbed  on  the  surface  and 
therefore  change  its  characteristics.  In  the  present  con¬ 
text  it  is  assumed  that  the  effective  surface  is  invariant 
with  incident  gas  properties.  For  two  derivations  of 
reciprocity,  both  classical  and  quantum  mechanical;  see  the 
paper  by  Kuscer  [26]. 

Equations  1,  2  and  4  are  the  basic  relations  that 
must  be  satisfied  by  the  scattering  kernel.  Notice  that 
Equations  2  and  4  are  linear  in  the  scattering  kernel. 

Thus,  if  P1  and  P2  satisfy  Equations  1,  2  and  4  then  P3  = 

0  Pj  +  (1-P)  P2  satisfies  Equations  2  and  4  for  all  values 
of  (3  and  satisfies  Equation  1  for  0  <  £3  <  1.  On  the  other 

3If  the  gas  adjacent  to  a  wall  were  in  such  a  state 
of  non-equilibrium  that  the  equilibrium  of  the  wall  were 
destroyed,  then  Equation  4  would  no  longer  be  valid,  e.g., 
for  extremely  high  rates  of  heat  transfer. 
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hand,  if  and/or  P2  do  not  satisfy  Equation  1  it  may  be 
possible  for  some  values  of  p  that  Pg  satisfies  Equations 
1,  2  and  4.  Therefore,  in  the  analysis  that  follows  in¬ 
dividual  scattering  kernels  are  not  disqualified  because 
of  failure  to  fulfill  Equation  1  since  a  linear  combination 
with  others  may. 

It  can  also  be  shown  that  if  P^  satisfies  Equation  4 
then  P2  also  satisfies  the  reciprocity  relation,  where  P2  is 
given  by: 

P2(v',v)  =  P1[2(T'-  J)  j-v’,^]  (5) 

Equation  5  is  simply  a  reflection  of  P2  about  the  normal 
velocity  and  it  is  easy  to  show  that  Equations  1  and  2 
are  also  fulfilled.  Thus,  if  a  forward  scattering  kernel 
satisfies  reciprocity  then  a  backscattering  kernel  also 
satisfies  reciprocity. 

1.3  SOME  SCATTERING  KERNELS  IN  VIEW  OF  RECIPROCITY 

1.3.1  Complete  Accommodation 

The  scattering  kernel  for  complete  accommodation  of 
a  gas  to  a  surface  is  the  same  as  the  effusive  flow  of  a 
gas  at  T  as  assumed  by  Maxwell  [3].  The  kernel  becomes 

V(v'  ,v)  =  (2/tt)  (2RTw)"2('v’-T)  exp(-v2/2RTw)  (6) 

1.3.2  Specular  Reflections 

The  scattering  kernel  for  specular  reflection  is: 
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(7) 


P(v*,v)  =  6(v  '-v  )6(v  *+v  )6(v  ’-v  ) 

\  /  v  x  x7  v  y  y'  v  z  z' 

where:  6(t)  is  the  Dirac  delta  function  and  has  the 

following  properties: 

6(t)  =  0 ,  t  ¥  0 

/  <5  (t)  =  1 


1.3.3  Diffuse  Elastic  Scattering 

The  scattering  kernel  for  diffuse  elastic  scatter¬ 
ing  from  a  surface,  also  called  Lambert  scattering  by 
analogy  with  light  reflected  from  a  rough  surface,  is: 

P( v"1  ,  v)  =  v  6(v'+v)/(ttv3)  (9) 

J 

1.3.4  Back  Scattered  Model 

Berman  and  Maegley  [27]  have  recently  (1972)  pro¬ 
posed  a  back  scattering  model  in  order  to  explain  the 
experimentally  observed  fact  that  the  flow  of  a  rarefied 
gas  through  a  tube  may  be  less  that  predicted  for  complete 
accommodations.  Many  years  earlier  Gaede  [28]  in  1913 
had  observed  the  same  phenomenon  and  concluded  that  based 
upon  Maxwell's  model  that  the  amount  of  specular  reflec¬ 
tivity  should  be  negative.  The  back  scattering  model  is: 

P(T' ,v)  =  6(v'+v  )6(v  '  +  v  )S(v  '+v  )  (10) 

vj/  vxx  yy  zz7 

This  expression  can  be  obtained  from  Equation  7 
by  the  application  of  Equation  5.  The  microscopic  observation 
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of  back  scattering  is  very  difficult  since  any  molecular 
detector  tends  to  disturb  the  incident  molecules. 

1.3.5  Epstein's  Scattering  Model  [22] 

As  mentioned  earlier,  Epstein  f 22 1  has  proposed  an 
empirical  scattering  kernel.  His  expression,  in  terms  of 
present  nomenclature,  is: 

PCT*  ,”v)  =  v  exp(-v2/2RT  )F(Tt )F(v)/[  /  v  exp(-v2/2RT  )• 
y  w  v  ^0  y 

y 

F(v) d3v]  +  6(v  ’-v  )6(v  ’+v  )6(v  ’-V  )  [l-F(T’)  ]  (11) 

ax  y  y  Z*  & 

By  substitution  it  can  be  shown  that  Equation  11  satisfies 
Equation  4  if  F(v)  =  F(-:v)  and  satisfies  Equation  1  if 
0  <  F(v)  <  1,  otherwise  F(v)  is  an  arbitrary  function  and 
satisfies  Equation  2  for  any  F(v) .  For  the  limits  of 
F(v)  s  0  and  F("v)  -  1,  Equation  11  reduces  to  specular  and 
diffuse  scattering,  respectively. 

1.3.6  Nocilla's  Scattering  Model  [19,  20] 

Nocilla’s  normalized  reflection  model  in  terms  of 
present  nomenclature  is: 

PCv-’  ,~v)  =  (2/7r)  (2RTr)-2  vy  exp[-(v-u).2/2RTr]/[(7r)1/2(uy/ 
^2RTr)  {l+erf  (uy/^/2RTr)}  +exp  (-uy2/2RTr>  ]  (12) 

where:  Tr  =  Tr  (‘v’* ) 

u  =  u(\T’ ) 
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The  quantities  Tr  and  u  are  the  temperature  and  mean  velocity 

of  the  assumed  reemitted  Maxwellian.  For  T  =  T  and  iT  =  0 

r  w 

Equation  12  reduces  to  Equation  6  and  for  Tr-»  0  and  iT  = 
v'-2(v‘‘j)T  Equation  12  approaches  Equation  7.  It  is  diffi¬ 
cult  to  ascertain  if  Nocilla's  model  satisfies  Equation  4. 
However,  Kuscer  [7]  has  stated  that,  in  general,  Equation  12 
does  not  satisfy  reciprocity. 

1.3.7  Summary 

By  direct  substitution  and  inspection  it  can  be  seen 
that  Equations  6,  7,  9  and  10,  as  well  as  Equation  11 
satisfy  the  reciprocity  condition  as  expressed  in  Equation  4 
as  well  as  satisfying  Equations  1  and  2.  Note  that  all 
these  kernels  except  Equation  6  are  discontinuous  in  velocity 
space  and  would  therefore  not  be  suitable  for  describing  a 
real  microscopic  gas-surface  interaction.  It  is  widely 
recognized  [24,  25,  26]  that  Equations  1,  2  and  4  are  neces¬ 
sary  conditions  which  must  be  satisfied  by  the  scattering 
kernel  for  any  surface. 

Any  physical  model  of  the  gas-surface  interaction 
such  as  those  given  in  [14,  15  and  16],  is  solved  exactly, 
would  satisfy  Equations  1,  2  and  4.  However,  when  approxi¬ 
mations  are  made,  as  is  usually  done,  the  results  do  not 
necessarily  fulfill  reciprocity,  although  Equations  1  and  2 
can  be  satisfied.  While  agreement  with  data  can  be  achieved 
over  a  limited  range  by  a  suitable  choice  of  model  parameters, 
these  models  cannot  be  extended  with  confidence  outside  the 
range  of  known  validity. 
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1.4  COORDINATE  SYSTEMS 

In  the  analysis  that  follows,  the  scattering  kernel 
is  solved  as  a  product  solution  in  three  different  coordinate 
systems;  rectangular,  cylindrical  and  spherical.  After  giving 
the  relations  among  the  various  coordinate  systems  the  three 
developments  will  be  presented  in  Section  II. 

Referring  to  Figure  1,  Appendix  A,  the  coordinates 
in  each  system  are: 


Coordinate  System 
rectangular 
cylindrical 
spherical 


Coordinates 


x,y,z 


p,  0,  cp 


The  cylindrical  and  spherical  systems  are  not  identical  with 
those  normally  encountered.  These  coordinates  were  chosen 
to  match  those  commonly  used  by  gas-surface  experimentalists 
[18,  29]. 

The  relations  between  velocity  components  and  angles 
in  these  coordinate  systems  are  shown  in  Table  I,  page  18. 
Also  given  is  the  expression  for  an  incremental  volume  in 
velocity  space. 
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TABLE  I 

RELATIONS  BETWEEN  VELOCITY  COMPONENTS 


Coordinate 

System 

Rectangular 

Cylindrical 

Spherical 

v  = 

X 

V 

X 

vrcos ^ 

VpCos<psin0 

Rectangular 

V  = 

y 

vy 

V 

y 

VpCOsqicosB 

vz  = 

vz 

vrsin^ 

v  sin© 

P 

V  = 

V 

V 

V  cos ©cos 0 

y 

y 

y 

p 

Cylindrical 

n 

> 

\l 

/v  2+v  2 
f  X  z 

vr 

Vp\Jl ”  C  OS  2  <pc  os  2  0 

if  = 

Tan”1 ( v  / v  ) 

Z  Zi 

if 

Tan”1  (Tanqp/sin0) 

V  = 

s 

2+ V  2+v  2 

L  2+v  2 

V 

p 

V  x  y  z 

v  r  y 

P 

Spherical 

0  = 

Tan-1 ( v  /v  ) 

a  y 

Tan-1  (  Vj.cosif'/Vy 

0 

<p  = 

Tan-1 

V  /  fv  2+v  2 

z  v  x  y 

o  •  -1 

Sin 

Jvrsin^/y/vr2+vy2J 

<P 

Incremental 

Q 

Volume  in 
Velocity  Space 

d  v  = 

dv  dv  dv 
x  y  z 

v  d^dv  dv 
r  r  y 

v  2cosqpdv  d0d^ 

P  p 
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II.  THEORETICAL  DEVELOPMENT 

Some  mathematical  solutions  to  Equations  1,  2  and  4 
will  be  obtained  in  this  section.  The  questions  may  arise 
as  to  how  a  scattering  kernel  can  be  determined  without  con¬ 
sidering  the  physics  of  the  gas-surface  interactions.  The 
situation  is  similar  to  that  in  a  gas  at  equilibrium,  de¬ 
tailed  balancing  is  valid  and  the  molecular  velocity  distribu 
tion  is  obtained  independently  of  the  physical  interaction 
between  molecules.  At  the  wall,  however,  a  molecule  inter¬ 
acts  for  only  a  finite  length  of  time  and  therefore  accom¬ 
modation  is  incomplete.  In  any  case,  the  set  of  functions 
that  satisfy  Equations  1,  2  and  4  contains  all  solutions  to 
any  physical  model  of  the  surface. 

Borrowing  from  Maxwell  the  idea  that  each  velocity 
component  in  a  gas  is  independent,  it  is  assumed  that  the 
interaction  of  each  velocity  component  with  the  surface  is 
independent.  This  is  equivalent  to  assuming  that  P(  aT*  ,~v) 
can  be  represented  by  a  product  of  independent  probability 
functions,  one  for  each  coordinate  direction. 

2.1  RECTANGULAR  COORDINATE  SYSTEM 

As  has  been  previously  stated,  it  is  assumed  that 
the  scattering  kernel  Pfv-’  ,"v)  can  be  written  as  a  product 
solution  in  rectangular  coordinates  as  follows: 

P(v',v)  =  Px(vx',vx)Py(vy’,vy)Pz(vz',  vz)  (13) 
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Substituting  this  expression  into  Equations  2  and  4 
and,  since  the  coordinates  are  independent,  equating  the 
terms  containing  each  coordinate,  gives  for  the  normaliza¬ 
tion  and  reciprocity  relations,  respectively: 

/  p ^ C vx *  i  vx)  dvx  =  1  (14-a) 


f  P  ( v  ’ ,  v  )  dv  =1 
J  yv  y  y'  y 

0 


00 


/  dvz  =  1 


(14-b) 


(14-c) 


and: 

exp(-vx,2/2RTw)Px(vx,J  vx)  =  exp(-vx2/2RTw)Px(-vx,-vx’) 

(15-a) 

-vy’exp(-vy,2/2RTw)Py(vy’,vy)  =  vy  exp(-vy2/2RTw> Py (-vy , -vy ’ ) 

(15-b) 

exp(-vz,2/2RTw)Pz(vz',vz)  =  exp  (- vz2/2RTw>  Pz  (-  vz, -vz  ' ) 

(15-c) 

Notice  that  the  equations  for  x  and  z  are  identical,  there¬ 
fore  only  the  analyses  for  x  and  y  will  be  presented  in 
detail. 

From  a  consideration  of  Equations  1  and  13  it  is 
obvious  that  P  (v  *,v  ),  P  (v  1  ,v  ),  and  P  (v  ’ ,v  )  must  be 

X  X  9  x  9  y  '  y  '  y7  *  Z  Z  9  Z^ 
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non-negative.  New  functions  are  defined  as  follows: 


f  (v  ' , v  )  =  Pv(v  ' , v  ) exp [ (A  2-l)v  ,2/2RT  +A  2v  2/2RT  ] 
X  X  X7  xvx  x  x  7  X  wxx  WJ 


(16-a) 


f  ( v  ' ,  v  ) 

y  y  y 


( V  i  V  )  09  9  9 

i - Z—  exp  [(A  -l)v  '  /2RT  +A  v  /2RT  ] 

vy  ^ 1  v  y  y  wyy  wJ 


(16-b) 


where  A  and  A  are  constants  equal  to  or  greater  than  unity, 
x  y 

f  and  f  were  defined  to  simplify  the  reciprocity  relation, 
x  y 

Substituting  Equation  16  into  Equations  14  and  15 
yields  for  the  normalization  conditions: 


!  fx(vx',vx)exp[-(Ax2-l)vx'2/2RVAx2vx2/2RTw]  dvx  >  1 


(17-a) 


/  V  f  (v  \  v  )  exp  [-  (A  2-l)v  '  2/2RT  -A  2v  2/2RT  ]  dv  =  1 
J  y  y  y  y  y  y  wyy  wJ  y 


(17-b) 


and  for  the  reciprocity  relation: 


f  (v  ' ,  v  )  =  f  (-v  ,-v  ’) 

X  X  X  X  X  X 


(18-a) 


f  ( v  ' , v  )  =  f  (-v  ,-v  ' ) 

y  y  y  y  y  y 


(18-b) 


From  Equations  1  and  16  it  can  be  seen  that  f  (v  ' ,v  )• 

A  A  A 

and  fy(vy,vy)  must  also  be  non-negative. 

Let  f  (v  1 ,v  )  and  f  (v  ’ ,v  )  be  represented  by 

xxx  yyy 

double  infinite  series  as  follows: 
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00 

00 

f  ( v  * ,  v  ) 

X  X  X 

■  Z 

j=0 

Z  alj  VX'‘  vxJ 

i=0 

(19-a) 

00 

CO 

f  ( V  *  ,  V  ) 

y  y  y 

-  z 

j=0 

V  b.  . v  ’ 1  v  ^ 

L.  y  y 

i=0 

(19-b) 

Reciprocity  requires  that 

a.  .  =  (-1) 
J-J 

i+ j  a  . . 
3  i 

(20-a) 

b.  .  =  (-1) 

i+j  b 

ji 

(20-b) 

Taking  advantage 

of  the  fact  that  reciprocity 

imposes 

a  certain 

kind  of 

symmetry  on  the  coefficients 

in  the 

expressions  for  fx  and 

f  ,  the  infinite  series  in  Equation 

Jr 

19  can 

be  rewritten  as  a 

double  sum  on  the  diagonals  of  the 

array  of  terms  as 

follows : 

00 

f  (v  \v  ) 

X  X  X 

-  z 

n=0 

An  Xn  (vx''vx> 

(21-a) 

00 

f  ( V  *  ,  V  ) 

y  y  y 

-  z 

n=0 

Bn  Yn  <vy''vy> 

(21-b) 

where: 

X  (v  ’,v  ) 
n  x  ’  x ' 

00 

-z 

k=0 

cn,k<vx'Vk[vx”+(-vx,)n1 

(22-a) 

22 
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W’V  "  I  dD,k{vy'vy>kfvyn+(-vy'>”]  (22*b> 


k=0 


with : 


_  ak>  k 

c0,k  257“ 


d  -  bk,k 
a0,k  2b^“ 


ak,  n+k 

n’k“-A 7“ 


d  t 

n,k  Bn 


^  n  =  1,  2, 


•  •  • 


By  inspection  it  can  be  seen  that  Equation  22  auto¬ 
matically  satisfies,  for  arbitrary  non-negative  integer 
value  of  n,  the  reciprocity  relation  as  expressed  in  Equation 
18.  The  only  restriction  upon  Equation  21  is  that  it  be 
non-negative  and  be  normalized  according  to  Equation  17. 

In  order  to  shorten  the  equations,  the  arguments  of 
variables  will  be  hereafter  omitted  unless  needed  for  clarity. 
Also,  velocities  will  be  non-dimensionalized  by  ^2RTw,  e. g. : 

v  - 

Sx  -  vx/y2RTw  (23) 

S  - 
etc. 
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It  is  now  assumed  that  each  and  Yn  is  individually 
normalized  by  Equation  17  as  follows: 


VO  dSx 


1  (24-a) 


(2RTW)  ”  SyYn  exp[-(Ay2-l)Sy'2  -Ay2Sy2]  dSy  =  1 

0  (24-b) 

Overall  normalization  is  satisfied  by  requiring  that: 


(25-a) 


00 

[b„-i  (25-*» 

n=0 

where  A„  and  B  are  arbitrary,  subject  to  the  above.  It  will 
n  n 

become  apparent  that  if  for  some  n,  Xn  or  Yn  is  not  normal¬ 
izable  then  the  corresponding  A  or  B  must  be  zero  in  order 

n  n 

to  preserve  the  overall  normalization. 

It  is  not  necessary  that  X  ,  Y  ,  A  or  B„  be  non- 

n  n  n  n 

negative  but  only  that  f  and  f  and  therefore  P  be  non- 

x  y 

negative. 

Each  X  or  Y  that  satisfies  the  above  conditions 
n  n 

leads  to  a  scattering  kernel  in  the  corresponding  coordinate 
direction  by  the  use  of  Equations  21  and  16.  Thus: 

Px,n  -  Xn  exp[-(Ax2-l)Sx’2  'AxVl  (26-a) 
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P  =  /2RT  S  Y  exp  [-(A  2-l)S  ,2  -A  2S  2]  (26-b) 

y ,  n  V  wyn  y  '  y  yyJ 


and: 


Xj  ^x,n 


n=0 


P  -  Y  B  P 
y  n  y,n 


n=0 


Letting; 


(27-a) 


(27-b) 


c 


n,  k 


Gn ,  k  V2»V#x> 


n+2k+l 


(28-a) 


and: 


dn,k  '  d„.k  ^“VV 


n+2k+2 


(28-b) 


and  substituting  in  Equation  22  gives: 


X„  = 


<2RT*)1/2  I  ax”+21!+1  Sn>kC Sx'sx)k  [sxn+<-sx')n] 


k=0 


(29-a) 


Y„  '  <2RVl  t  Ay"+2k+2  dn,k<Sy,Sy>k  t Syn+ <-Sy ’ > " ] 


k=0 


(29-b) 


2.1.1  The  Evaluation  of  Px  >n 

Substituting  Equation  29-a  into  Equation  24-a  yields: 
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oo 

oo 


1  I  V+2k+1  vk<VVk[s*"+<-V>nJ  "•w.V 


dSx  - 


k=0 


exp[(Ax2-l)Sx’2] 


(30) 


Integrating  term  by  term  on  the  left  hand  side  where  each 
term  is  of  the  form: 


00 


/  tra  exp(-p2t2)  dt  = 


0  if  m  is  odd 


m+1 


f(m/2  +  1/2) /p  if  n  is  even 


(31) 


and  expanding  the  right  hand  side  in  an  infinite  series, 
gives  the  following  power  series  in  S’: 


V  A  2k  £  S  ,2k 
L  x  c2n,2k  ^x 

k=0 


r(k  +  n  + 


i/2)  -  £ 

k=0 


2n+2k  ~ 
x  c2n , 2k* 


oo  (A  2_^k  „  ,2k 

(S  ’) 2k+2n  r(k  +  1/2)  =  y  __x - - x -  (32) 

X  i  K  ■ 

k=0 

for  even  indices  of  c„  ,  .  Referring  to  Equation  30,  for  the 
first  index  and  the  second  index  odd  the  coefficients  are 
arbitrary  since  the  integrals  and  the  right  hand  side  of  Equa¬ 
tion  32  are  zero.  If  the  first  index  is  odd,  the  integrated 
form  of  Equation  30  cannot  be  satisfied,  because  there  are 
only  odd  powers  of  S  ’  on  the  left  hand  side  and  only  even 

X 

powers  on  the  right  side.  Therefore,  ^2n+1(vx' , vx)  cannot  be 
normalized  and  A2n+1  must  be  zero.  The  only  case  that  must 
be  investigated  is  that  given  by  Equation  32. 
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Equating  the  coefficients  of  like  powers  of  Sx'  in 
Equation  32  yields: 


cn ,  k 


_ 1 _ 

(k/2)  !  r(k/2  +  n/2  +  1/2) 


n 

.  ,  /Ax  ~1\2  _ 1 _ 

n,k  \  Ax2  I  <k/2>;  r(k/2  +  n/2  +  1/2) 

f  r(k/2  -  n/2  +  1/2)  k  n 

n,k-n  r (k/2  +  n/2  +  1/2) 


k  <  n 


(33) 


for  k,n  even.  For  odd  k,  c  .  is  arbitrary.  However,  it  is 

n  f  k 

assumed  that  Equation  33  also  gives  the  odd  coefficients. 
Given  values  of  A  and  even  n,  and  using  Equations  26-a, 

X 

29-a  and  33,  P  can  be  determined  as  a  function  of  S' 

'  x,  n  x 

and  S  .  The  case  of  n=0  can  be  solved  analytically.  For 

X 

n=0  Equation  33  becomes: 


1 

2 


_ 1 _ 

r(k/2  +  1)  T(k/2  +  1/2) 


(34) 


using  the  fact  that  (k/2) J  =  r(k/2  +  1).  Equation  6.1.17 
of  [30]  permits  writing: 

T(k/2  +  1)  T(k/2  +  1/2)  =  —  +1J-  =  ^  k:  (35) 

2k  2k 


Substituting  Equation  35  in  Equation  34,  Equation 
29-a  becomes: 
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,/9  (A  -1)  A  . 

XQ  =  (27TRTw)"i/2  Ax  ^  ^ - (2SX»  Sx)K  (36) 

k=0 

By  inspection  it  can  be  seen  that  this  is  the  expan¬ 
sion  of  the  exponential  function,  therefore  Equation  36 
becomes; 

X„  =  (2ttRT  )“1/2  a  exp [2 (A  2-l)1/2  A  S  *S  ]  (37) 

O  W  X^lvX  XX  XJ 

Substituting  in  Equation  26-2  gives: 

P*,0  -  ^'1/2  Ax  exp  -  {(Ax2-1>1/2  V  -Aa)  2  <38> 

This  equation  contains  most  of  the  experimentally  observed 
properties  of  gas-surface  interactions.  For  the  case  of 
Ax=l  Equation  38  reduces  to  the  fully  accommodated  diffuse 
reflection  and  for  A  <»  approaches  6(S  '-S  )  which  is 
specular  reflection.  This  specular  limit  can  be  shown 
from  a  consideration  of  Equation  21.9-18b  in  [31].  It  will 
be  shown  later  (Section  3.3)  that  Ax  is  related  to  the 
thermal  accommodation  coefficient,  a,  and  is  in  fact,  for 
n=0,  given  by  the  following: 

Ax  =  (ax)1/2  (39) 

The  x  subscript  was  added  to  a  to  allow  the  accom¬ 
modation  coefficient  to  vary  with  coordinate  direction. 

Figure  2,  Appendix  A,  illustrates  the  results  of  Equation  38 
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with  /2RT  P  n  given  as  a  function  of  S  ,S  1  and  a  . 
y  w  x,  u  xx  x 

For  n  greater  than  zero  (n=2, 4, 6, . . . )  it  was 
necessary  to  solve  Equations  26-a,  29-a  and  33  numerically. 

A  computer  program  was  written  to  solve  these  equations. 

The  results  are  shown  for  J2RT  P  in  Figures  3  and  4, 

V  W  X  |  II 

Appendix  A,  for  n=2  and  4,  respectively. 

2.1.2  The  Evaluation  of  Pyin 

Substituting  Equation  29-b  into  Equation  24-b  yields: 


oo  r—> 

L  sy  L 


Ay"+2k+2  Jn,k(S/Sy)lt[Syn-(-Sy')nlexp(-Ay2sy2>  dSy  “ 


k=0 


exp[(Ay2-l)Sy’2]  (40) 

As  was  done  for  the  x  coordinate  the  left  hand  side 
can  be  integrated  term  by  term.  Each  term  is  of  the  form: 


/  tm  exp(-p2t2)  dt  =  r(m+l/2) /2pm+1  (41) 

0 


Performing  the  integration  and  expressing  the  right  hand 
side  as  an  infinite  series,  Equation  40  becomes: 

\  E  <-»"  Ay”+k  Sn,k  Sy’k+D  r<k/2  +  l)  + 
k=0 

5  fAyk  3n,k  sy'k  rCk/2  +  n/2  +  »  ' 

k=0 
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oo 


I 


(A  2-l)k  S  ,2k 
y  y  y 

k.f 


(42) 


Since  this  equation  holds  for  arbitrary  value  of  Sy’  the 
coefficient  of  like  power  of  can  be  equated,  and  solving 
for  dn  k  yields: 


dn,  k 


_ 2 _ 

r(k/2  +  n/2  +  l)(k/2).' 


0 


k  even,  <  n 


k  odd,  <  n 


_ 2 _ 

r(k/2  +  n/2  +  l)(k/2).' 


(-1)" 


d  , 
n ,  k-  n 


P(k/2  -  n/2  +  1) 
r(k/2  +  n/2  +  1) 


k  even,  >  n 


n 

d  ,  =  -  (-1)  d  . 
n,k  v  n,k-n 


T(k/2  -  n/2  +  1), 
T(k/2  +  n/2  +  1) 


odd,  >  n 


(43) 


Using  Equations  26-b,  29-b  and  43  Py  n  2RTw  can  be 

determined  as  a  function  of  S  ,S  ’  and  A  for  given  values 

y  y  y 

of  n.  As  was  true  for  P  n,  P  n  can  be  obtained  in  closed 

x , u  y,u 

form.  For  n=0  Equation  43  becomes: 


r(k/2  +  1)  (k/2) 


for  k  =  0,  2,  4,  .  .  . 


dg  k  =  0  for  k  "  1,  3,  5,  ...  (44) 
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Substituting  Equation  44  in  Equation  29-b  yields, 


00 


Y0  -  CRT  )-1  l  <Ay2-l,k  Ay2* 

kt'o  y  y  r(k+i)k: 


2k 


(45) 


From  inspection  (see  Equation  9.6.10  in  [30])  it  can  be  seen 
that  Equation  45  is  the  expansion  of  the  modified  Bessel  func¬ 
tion  of  the  first  kind,  I  .  Equation  45  therefore  becomes: 

=  A  2  (RT  )-1  In  (-2 /a  2-l  A  S  'S  )  (46) 

o  y  w  ovvy  yyy 

Substituting  Equation  46  in  Equation  26-b  gives:4 


Sy’V 

exp[-(Ay2-l)  Sy ’ 2  -Ay2Sy2]  (47) 


As  was  the  case  for  P  n,  Equation  47  contains  most 
of  the  experimentally  observed  properties  of  gas-surface 
interaction.  For  the  limit  of  Ay=l  Equation  47  reduces  to 
the  fully  accommodated  diffuse  reflection  and  for  Ay 
using  the  asymptotic  expansion  for  Iq,  Equation  47  reduces  to: 


6Xp[-Ay2(Sy '  +  Sy)2] 


(48) 


or  as  was  shown  for  P  n 

X  )  vs 


4After  the  present  analysis  was  complete,  it  was 
found  that  Kuscer  et  al  [7]  had  obtained,  in  a  completely 
different  manner,  a  result  similar  to  that  given  in 
Equations  38  and  47. 
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(49) 


Py,0  V2RTw  ~6<Sy'  +  V 


Thus  for  A  —•»«>,  S  =  -S  * ,  which  is  also  specular  reflec- 

y  y  y 


tion . 


Figure  5,  Appendix  A,  illustrates  the  results  of 


Equation  47  with  /2RT  P  0  given  as  a  function  of  S  ,  S  ' 

v  w  jr  >  u  y  y 

and  ay.  For  n=l,  2,  Equations  26-b,  29-b  and  43  were  solved 

numerically.  The  results  for  P„  are  shown  in  Figures  6 

y ,  n 

and  7,  Appendix  A,  for  n=l  and  2,  respectively,  as  a  function 
of  S  ,  S'  and  a  where: 


a  =  A 

y  y 


(50) 


2.2  CYLINDRICAL  COORDINATE  SYSTEM 


Keeping  in  mind  the  cylindrical  coordinates  y,  r, 
and  i/  as  given  in  Figure  1,  Appendix  A  and  Table  I,  page  18, 
the  scattering  kernel  is  now  assumed  to  be  representable 
as  follows: 


P(v\v)  =  Py ( vy  1  > vy )  pr<V>V  (51) 


In  cylindrical  coordinates  Equation  2  becomes: 


7T  oo  oo 

III  pfv’1  *v)  v  dv  dv  d^  -  1 

-■n  o  o  r  y  r 


(52) 


Substituting  Equation  51  into  Equation  52  the 
following  are  obtained  for  the  normalization  relation  for 
each  coordinate,  respectively: 
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/  P„  dv  =  1  (53-a) 

0  y  y 

7  vr  Pr  dvr  =  1  (53-b) 

0 

7 r 

/  P ,  dtf  =  1  (53-c) 

-7T  * 

The  reciprocity  relation  given  by  Equation  4  becomes  in 
cylindrical  coordinates  with  the  applications  of  Equation  51: 

-v  '  exp (—  S  ,2)  P  (v  *,v  )  =  v  exp (- S  2)  P  (- v  ,-v  ') 

y  y  y  y  y  y  y  y  y  y 

(54-a) 

exp(-Sr'2)  Pr(vr’,vr)  -  exp(-Sr2)  Pr(-vr,-vr')  (54-b) 

Px^i’^  =  (54-c) 

Notice  that  Equations  53-a  and  54-a  are  identical  to 
Equations  14-b  and  15-b,  respectively.  Therefore,  the 
results  for  the  y  coordinate  in  cylindrical  and  rectangular 
coordinates  are  the  same. 

2.2.1  The  Evaluation  of  Pr>n 

A  function,  f„(v  ',v„),  is  defined  as  follows: 
i*  r  r 

fr(vr’.vr)  =  Pr  exp[(Ar2-l)Sr'2  +  Ar2  Sr2]  (55) 

Utilizing  this  relation,  Equations  53-b  and  54-b  become, 
as  a  function  of  fr,  respectively: 
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/  vrfr  exp[-(Ar2-l)Sr’2  -Ar2  sr2]  dvr  -  1  (56) 

fr(vr’,vr)  =  fr("vr,"fr')  (57) 

These  two  equations  are  the  same  as  Equation  17-b  and  18-b 
except  for  the  subscripts.  This  leads  to  the  conclusion;: 

VV'V  $  fy^vr'  ’  vr)  (58) 

Assuming  that  Pr  is  given  by  a  linear  combination  of  in¬ 
dividual  kernels  as  was  done  in  Equation  27,  permits 
writing: 


00 


=  y  c  p 

n  r,  n 


» (59) 


n=0 


Subject  to,* 


E 


(60) 


n=0 


where,  from  a  consideration  of  Equations  58,  55  and  16-b; 


Py  n(vr,'vr) 

P  (v  * , v  )  =  — — - 

r,nv  r  ’  r'  v_ 


(61) 


For  n=0,  in  view  of  Equation  47,  Equation  61  becomes: 


(2RTw)Pr,0=2Ar2  Io(2Ar'/r2-1  W  exP[-<Ar2-1>Sr'2-Ar2Sr2] 

(62) 
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The  results  of  this  equation  are  shown  in  Figure  8,  Appendix 
A.  The  numerical  calculation  for  n=l,  and  2  are  shown  in 
Figures  9  and  10,  Appendix  A,  respectively. 

2.2.2  The  Evaluation  of  P^, 

From  a  consideration  of  Equation  62  and  Figures  8, 

9  ,  and  10,  Appendix  A,  it  can  be  seen  that,  in  general, 

P  is  finite  at  v  =0.  In  view  of  this  and  since  for  normal 
r ,  n  r 

reflected  velocities  (vr=0)  the  reflected  probability  is 
independent  of  the  only  solution  is  for  to  be  equal  to 
a  constant.  From  Equation  53-c  it  is  found  that: 


which  also  satisfies  Equation  45-c. 

2.3  SPHERICAL  COORDINATE  SYSTEM 

Using  the  spherical  coordinates  as  given  in  Figure  1, 
Appendix  A,  and  Table  I,  page  18,  the  scattering  kernel  is 
now  assumed  to  be  representable  as  follows: 

P  ( v  * ,  v)  =  P  (v',v)  P^CS^S)  (64) 

Note  that  the  p  subscript  is  not  needed  on  the  velocity  com¬ 
ponents  since  the  velocity  vectors  coincide  with  the  coor¬ 
dinate  direction  and  the  velocity  component  is  therefore 
equal  to  the  molecular  speed  when  the  velocity  is  directed 
away  from  the  surface  and  to  its  negative  when  directed 


35 


AEDC-TR-72-1 58 


toward  the  surface.  In  spherical  coordinates  Equation  2 
becomes : 


so  7r/2  7T/2 


/  /  / 

0  — 7r/2  -7T/2 


P  v  cos  <p  dq)  d0  dv  =  1 


(65) 


Substituting  Equation  64  into  Equation  65  the  following 
normalization  relations  are  obtained  for  each  coordinate, 
respectively: 

x  o 

f  P  v  dv  =  1  (66-a) 

0  P 


7T/2 

/  P9dS 
-tt/2 


1 


7r/2 

/ 

-7r/2 


P^  cos  cp  dcp  =  1 


(66-b) 


(66-c) 


In  spherical  coordinates,  after  the  application  of 
Equation  64,  Equation  4,  the  reciprocity  relation,  becomes: 

-v'  exp(-S’2)  P  (v',v)  =  v  exp(-S2)  P  (-v,-vf)  (67-a) 

P  P 

cos  ei  Pd(eite)  =  cos  e  Pe(fl,P.)  (67-b) 

cos  qpi  ,  <p)  =  cos  cp  P^c cp,  cpi)  (67-c) 
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From  Equations  1  and  64  it  can  be  seen  that  P  ,  P_ 

P  y 

and  P^  must  be  non-negative.  This  requirement  along  with 
Equations  66  and  67  are  the  conditions  that  must  be  satis¬ 
fied  by  P  ,  Pa  and  P  . 

p  9  9 

2.3.1  The  Evaluation  of  P p  n 

The  development  of  the  solution  for  P^  will  closely 

parallel  the  solution  for  P  and  P  .  A  function  f  is 

x  y  p 

defined  as  follows: 


fp(v’’v)  =  /  exPt  (Ap^-l)  S’2  +  ApZ  SZ] 


2  c.2 


(68) 


Equations  66-a  and  67-a  become,  in  terms  of  f^,  respectively: 


f  v3  f  exp[- (A  2-l)  S'2  -A  2  S2]  dv  -  1 
0  P  P  P 


(69) 


and : 


f  (v'.v)  =  f  (-v,-v() 

r  F 


(70) 


It  is  now  assumed  that  f  can  be  represented  as: 


f  =  y  D  w_ ( v' , v) 

p  i—t i  n  n 


n=0 


(71) 


where : 


Wn ( v '  ’  v)  =  Z  gn,k  <v’v)k[vn  +  (-v’)n]  (72) 

k=0 
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and  the  Dn's  are  arbitrary  constants  subject  to: 


00 


n=0 

Notice  that  the  assumed  form  for  and  therefore 

n 

f  as  expressed  in  Equations  71  and  72  automatically  satis- 

H 

fies  the  reciprocity  relation  given  in  Equation  70.  There¬ 
fore  the  only  constraints  upon  f  are  the  normalization 

P 

conditions  given  by  Equation  69  and  the  condition  of  being 
non-negative. 

Assuming  that  is  given  by  a  linear  combination 
of  individual  kernels  as  was  done  in  Equation  27,  yields: 


>  =  y  d  p 

p  n  p,n 

n=0 


(74) 


Equation  72  becomes,  in  terms  of  S  and  S' 


\  -  I  *n,k  <S'S)k[Sn  ♦  <-S')°]A 


2k+n+4 


k=0 


(75) 


where: 


y2 


RT 


w 


n,k  &n, k  \  A, 


n+2k+4 


(76) 


For  Dn  =  1  and  =  0  (j  ^  n)  Equations  71,  74  and 
75  may  be  used  along  with  the  dimensionless  speeds  to  permit 
writing  Equation  68  as  follows: 
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(2RTw)3/2  Pp  n  =  S  exp[-  (Ap2-1)  S 1  2  -Ap2  S2  ]  * 


Z  ®n,k 
k=0 


(S’S)k[Sn  + 


(-S’)n] 


A 


n+2k+4 

P 


(77) 


Substituting  Equations  71  and  72  into  Equation  69  gives: 


00  q  r — ' 

n+2k+4  f  S  )  ~  ,.,..kIon  ,  -..n^  ,  .  2  _2S 

i  gn>k(s  S)  [S  +  (-S’)  J  exp(-A^  S  )  dS 


k=0 


=  exp[  (Ap2-1)  S'2] 


(78) 


Equation  78  can  be  integrated  term  by  term  on  the  left  hand 
side  where  each  term  is  of  the  form  given  by  Equation  41. 
Performing  this  integration  and  expanding  the  right  hand 
side  as  an  infinite  series,  Equation  78  becomes: 


|  I  Apk  2n,k  Vk  r(k/2  +  D/2  +  2)  + 
k=0 


(-1)° 

2 


Z  Apk+m  ^n,k  (S’>k+n  r(k/2  +  2)  = 

k=0 


Y  (A72-l)k  S'2k  (79) 

k=0 

Since  this  equation  is  valid  for  arbitrary  values  of  S',  the 

coefficients  of  like  powers  of  S'  can  be  equated.  Solving 

this  set  of  equations  for  g  .  gives: 

n  i  k 
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gn,k 


_ 2 _ 

(k/2) :  r(k/2  +  n/2  +  2) 


k  even,  <  n 


gn,  k 


0 


k  odd,  <  n 


gn,k 


_ 2 _ 

(k/2):  r(k/2  +  n/2  +  2) 


(-I)” 


gn, k-n 


r(k/2  -  n/2  +  2) 
r(k/2  +  n/2  +  2) 


k  even,  >  n 


gn,  k 


-(-Dn 


gn , k-n 


r(k/2  -  n/2  +  2) 
r(k/2  +  n/2  +  2) 


k  odd,  >  n  (80) 


Using  Equations  77  and  80,  P  (2RT  ) 3//2  Can  be 

p  i  n  w 

determined  as  a  function  of  S,S’  and  A  for  given  values  of 

P 

n.  As  was  the  result  for  Pv  n,  P_  _,  P  can  be  obtained 

x,  u  r,U  p , U 

in  closed  form.  For  n=0  Equation  80  becomes: 


g0,  k 


_ 1 _ 

(k/2).’  r(k/2  +  2) 


even  k 


g0,  k 


=  0 


odd  k 


(81) 


Substituting  Equation  81  into  Equation  77  yields: 

<2RTw)3/2  Pp0  -  S  exp[-(Ap2-l)S'2  -Ap2  S2]- 
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(S’S)2k  A  2k 


4  V"1  2  k  An 

2A  V  (Ad  _1)  - 9~ 

9  9  k.’  r(k+2) 


(82) 


From  Equation  9.6.10  of  [30]  it  can  be  seen  that 
this  last  sum  is  equal  to: 


Il(-2S'S  A 


-s's  Vap  1 


Equation  82  therefore  becomes: 


(2RT  )3/2  P  n 

w  p ,  0 


exp  [-(A  2-l)S'2  -A  2  S2]* 
P  P 


!A P  -1  <-s*> 


I1(-2S’SAp^Ap2-l) 


(83) 


Equation  83  is  plotted  in  Figure  11,  Appendix  A, 
3/2 

with  (2RT  )  P  n  given  as  a  function  of  S  for  various 

W  Pi  « 

values  of  S'  and  a  ,  where: 

P 


a  -  — 5 
P  A  ^ 
P 


(84) 


For  n  greater  than  zero  (n=l,2)  Equations  77  and  80 

3/2 

were  solved  numerically.  The  results  for  (2RT  )  P  .  and 

w  p  i  A 

3/2 

(2RT  )  P  r»  are  shown  in  Figures  12  and  13,  Appendix  A, 
w  p,2 

respectively. 
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2.3.2  The  Evaluation  of  Pq 

It  is  assumed  that  Pg  is  given  by  the  following 
infinite  series: 

00 

Pg  =  cos  0  e^  cos  ‘[  j  C @ ^  +  0)]  (85) 

j=0 

This  form  was  chosen  because  each  term  on  the  right 
hand  side  automatically  satisfied  Equation  67-b.  Substi¬ 
tuting  Equation  85  into  Equation  66-b  gives  for  the  normal¬ 
ization  condition: 


tt/2 

/ 

-7r/2 


CO 


f  cos  0  ^  ej  cos  +  0) ]  d0  =  1 

-nr/O 


(86) 


Each  term  in  this  equation  is  of  the  form: 


7T/2 

J  cos  8  cos  [3(0^  +  0)]  d0 

-7T/2 


By  using  the  trigonometric  identity  for  the  cosine  of  the 
sum  of  two  angles  this  integral  becomes: 

7T/2 


J  cos  0[cos  j0^  cos  j0  -  sin  sin  j0]  d0 


-tt/2 


which,  after  integrating,  yields 


cos 


x  fsin(.j-l)7r/2  .  sin (.i+l)7r/2l 

(jei}  l  J-l  j+1  J 
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which  may  be  reduced  to: 

-2  cosCjd^)  COS  (jTT/2) 

3^1 


This  expression  is  zero  for  all  odd  values  of  j 
except  j=l  where  it  is  of  the  form  zero  over  zero  which 
must  be  evaluated  by  l'Hopital's  rule.  The  final  result  is: 


7T/2 


I 

-7T/2 


cos  0  cos  [  j  ( +  0)]  d‘0  =  2 


»  3  = 


0 


=  7T/2  COS  0 i 

-2(-l)J/2  cos  (j0.) 
(j2-l) 


,  j  =  1 

,  3  even 


=  0 


,  j  odd,  >  1 
(87) 


Since  the  integral  is  zero  for  j  odd,  >  1  the  corresponding 

e.  is  arbitrary  since  that  term  does  not  contribute  to  the 
J 

sum  in  Equation  86. 

Substituting  Equation  87  into  Equation  86  gives  the 
following: 


00 


2  e^  +  7r/2  e^  cos  0^  -  2 


e 


2k 


k=l 


(-1;  cos  2k0^ 
(4k2-l) 


=  1 


(88) 


By  inspection  it  can  be  seen  that  this  equation  is 
equivalent  to  the  equation  for  the  Fourier  expansions  of 
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l-(7r/2)  e,  cos  0.  over  the  period  -tt/2  to  7t/2.  The  Fourier 
series  is: 


00 

1  -tt/2  e^^  cos  =  aQ/2  +  ^  cos  2k  0j^  (89) 

k=  1 


where  the  Fourier  coefficients  are  : 


7T/2 

=  2/rr  J  (l-7r/2  e^^  cos  0^  cos  2k  ©i  d©i  (90) 
-tt/2 


Notice  that  the  terms  containing  sin  2k0i  are  missing  from 
Equation  89  because  the  left  hand  side  is  an  even  function 
of  0i. 

Equation  90  can  be  integrated  to  give; 


a 


a 


2k- 1 


a 


2k 


2(1^) 

0 

(-1)  k2ex 
(4k2-l) 


1,2,3, 


(91) 


Substituting  Equation  89  into  Equation  88,  using 
Equation  91  and  equating  like  coefficients  of  cos  2k0^  gives: 


e0 


(92) 


e2k  =  "el  ’  k  ^  1 


It  has  already  been  shown  that  the  odd  e's  greater  than  1 
are  arbitrary. 
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Utilizing  Equation  92.  Equation  85  becomes: 


P0  =  cos 


[1  CO 

^  +  e^^  cos  (7^+0)  -e.^  ^ 

k=l 


cos  2k(0^+0)+ 


L  ®2k+i  o°sC2k-i)(e1+0)l 

k=l 


(93) 


By  inspection,  it  can  be  seen  that  the  first  sum 
in  Equation  93  will  not  converge  by  the  normal  definition 
since  the  coefficients  are  constant.  It  can  be  seen  that 
the  value  of  the  sum  at  0^+0=O  is  unbounded  since  all  terms 
are  "in  phase"  at  this  point.  However,  at  other  values  of 
0^+  0 1  the  sum  oscillates  from  positive  to  negative  with  un¬ 
limited  frequency  and  over  any  interval  is  of  unbounded 
variation.  This  sum  is  indeed  a  badly  behaved  function. 
However,  in  the  mean,  this  sum  acts  similar  to  a  Dirac 
delta  function.  The  Fourier  coefficients  for  the  delta 
functions  will  now  be  determined. 

Any  real  function,  f (t) ,  defined  over  the  finite 
interval  -t/2  <  t  <  t/2,  defines  a  Fourier  series  if: 

t/2 

[f(t)|  dt  (94) 

-t/2 

is  bounded  (see  [31],  Sections  4. 11- 2b).  The  Fourier  series 
is  defined  as: 


45 


AEDC-TR-72-158 


00 

^  aQ  +  ^  ak  cos  (27tkt/T>  (95) 

k=l 


if  f(t)  is  symmetric  about  t=0,  the  expression  for  afe  is: 

t/2 

afc  =  2/t  J  f(t)  cos  (27rkt/T)  dt  (96) 

-t/2 

The  delta  function  S(t)  satisfied  the  above  condition.  In 
fact,  if  we  let  the  interval  be  -7r/2  <  t  <  tt/2,  expression 
94  becomes : 

7r/2 

f  |s(t)|  dt  =  1  (97) 

-t/2 

Equation  96  gives  for  the  Fourier  coefficients  of  the  delta 
functi on: 


7r/2 

ak  =  2//?r  J  <5(t)  cos  (2kt)  dt  (98) 

— tt/2 


-  2/t 

Thus,  the  Fourier  coefficients  of  the  even  cosine 
terms  are  constants;  this  is  true  also  for  the  coefficients  of 
the  first  sum  in  Equation  93.  Notice  that  the  Fourier 
series  given  by  Equations  95  and  98  does  not  converge  uni¬ 
formly  to  the  delta  function  over  any  subinterval. 
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Convergence  is  guaranteed  only  for  a  function  of  bounded  \ 
variation  (see  [31],  Section  4.11-4)  and  the  delta  function 
is  unbounded. 

While  it  is  true  that  the  Fourier  series  of  the  delta 
function  is  not  equal  to  the  delta  function  itself,  it  may 
be  that  they  both  satisfy  the  same  integral  relation  as 
given  in  Equation  86.  Therefore,  substituting  for  the  series 
of  even  terms  in  Equation  93  the  delta  function  and  re¬ 
writing  the  series  of  odd  terms  gives: 

ri+ei 

=  cos  0  -  +  e,  cos  (0.+0)  -  e-,  7t/2  6(0. +0)  +  G(0.+0) 

8  L21  1  1  1  1 

(99) 

where  G(0j  +  0)  is  the  sum  of  odd  terms  in  Equation  93  and  is 
an  arbitrary  function  subject  to  the  following  conditions  : 

G(  0)  =  G(-0) 

G(0)  =  -G(tt-0) 

7T/2 

cos  0  G( ©i+0)  d0  -  0  (100) 

-tt/2 

Equation  99  does  indeed  satisfy  Equation  66-b  for  arbitrary 
ej^  and  G  subject  to  Equation  100.  However,  ej^  and  G  are 
still  limited  such  that  Pg  is  non-negative.  It  will  be 
shown  in  Section  2.4.3  that  e^  and  G  are  further  restricted. 

2.3.3  The  Evaluation  of 

Similar  to  what  was  done  for  Pg  it  is  assumed  that 
P  is  given  by  the  following  series: 
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COS 


COS 


[^(cp^cp)  ] 


(101) 


The  reciprocity  relation  as  expressed  in  Equation  67-c  is 
satisfied  by  Equation  101  for  arbitrary  h.'s.  Substituting 
Equation  101  into  Equation  66-c  gives  for  the  normalization 
relation: 


tt/2 


/ 

-Tr/2 


cos 


h . 
J 


cos  [jCcpj  +  qO  ]  dcp  =  1 


(102) 


Each  term  in  the  infinite  series  contains  an  integral  of  the 
form: 


it/2 


/ 

-tt/2 


cos 


cp  cos  [  j  (cpi+qj)  3  dcp 


which  becomes,  after  using  the  trigonometric  identity  for 
the  cosine  of  the  sum  of  angles: 

7r/2  tt/2 


cos  j  cp 


cos  cp  cos  j cp  dp  -  sin  j 


j”i  j 


cos  cp  sin  jcp  dp 


-7T/2 


-7T/2 


The  second  integral  is  zero  from  symmetry.  Integrating  the 
remaining  term  gives: 

4  cos  j  Pj^  sin  (j  7t/2) 
j (4-j2) 
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This  expression  is  zero  for  all  even  values  of  j  except 
j=l,  2  where  it  is  of  the  form  zero  over  zero  in  which 
case  it  must  be  evaluated  by  l'Hopital's  rule.  The  final 
result  for  the  integral  is: 


7T/2 


f  cos2  cp  COS  [jCaK+tp)]  dcp 


-7r/2 


tt/2  , 

4  cos  qp^ 


=  7t/4  cos  2cpi 
=  4/15  cos  3cp^ 


j  =  0 

j  =  1 
j  =  2 
j  =  3 


=  0 


even  j,  >  4 


f.»<J+l)/2  4 


COS 


jcPL  > 


odd  j 
(103) 


Since  the  integral  in  Equation  102  is  zero  for 
even  j’s  greater  than  two,  the  corresponding  h  is  arbi- 

J 

trary  since  it  does  not  contribute  to  the  sum  in  Equation 

102. 

Substituting  Equation  103  into  Equation  102  for 
each  term  in  the  series,  gives: 


00 

77/2  hQ  +71/4  hg  cos  2cp^  -  ^  (-l)k* 

k=0 


4  cos(2k+l)qp^, 
2k+1  (2k+l) (2k-l) (2k+3) 


(104) 
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Equation  104  is  equivalent  to  the  Fourier  expansion  of 
l-rr/2  hQ  -7t/4  h 2  cos  2<p^  over  the  range  — 7r/2  <  <  7r/2. 

Because  there  are  only  odd  multiples  of  cp^  in  the  series 
in  Equation  104  the  Fourier  series  must  be  skew-symmetric 
about  ±tt/2.  The  Fourier  series  is: 


00 

l-ir/2  hQ  — 7r/ 4  h2  cos  2^  =  aQ/2  +  ^  ak  cos  kq^ 

k=l 


(105) 


where : 


ak- 


4/t r-  2h. 


k  h, 
k2-4 


odd  k 


=  0 


even  k. 


Equating  Equations  104  and  105  and  solving  for  h^  gives 
for  odd  k: 


hk  =  +4/ 7T  -  2hQ  -  k2  [1/tt  -  hQ/2  +  h2/4]  (106) 

A  necessary,  though  not  sufficient,  condition  for  the  con¬ 
vergence  of  P  ,  as  given  by  Equation  101,  is  for  h,  to  be 

(D  K 

bounded  as  k-^».  This  fact  requires  that  the  coefficient 
2 

of  k  in  Equation  106  be  zero,  i.e.: 

h2  =  2h0  "  4/?r  (107) 


Equation  106  therefore  becomes  : 
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hk  =  4/7T  “  2h0  »  even  k  (108) 

Utilizing  Equations  107  and  108  and  the  fact  that  the  even 
h^'s  for  k  greater  than  2  are  arbitrary,  Equation  101 
becomes : 

P^  =  cos  cp[hQ  -  (4/7r-2hQ)  cos  2(cpi+cp)  + 

00 

(4/7r-2hQ)  ^  cos  (2k-l)  (q^+qj)  + 

k=l 


00 

Yi  h2k+2  cos  (2k+2>  ('Pi+tP>I  (109) 

k=  1 

Notice  that  the  first  infinite  sum  in  Equation  109 
has  constant  coefficients  as  was  the  case  in  Equation  93. 
However,  the  series  in  Equation  93  contained  even  multiples 
of  the  angles,  while  Equation  109  contains  only  odd  multi¬ 
ples. 

In  Section  2.3.2  the  Fourier  coefficients  of  the 
delta  function  were  shown  to  be  equal  to  a  constant  for 
the  even  cosine  terms  and  zero  for  the  other  terms.  It  will 
now  be  shown  that  a  Fourier  series  for  the  delta  function 
can  be  developed  that  contains  only  odd  cosine  terms.  Let 
f(t)  be  defined  as  follows: 

f(t)  =  6(t)  -  6 ( t+7r)  -  6 ( t— tt)  (110) 
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By  inspection  it  can  be  seen  that  Equation  110  is  equivalent 
to  6(t)  over  the  range  -77  <  t  <  7r.  From  Equation  96  with 
t  =  27 r  and  using  Equation  110  there  is  obtained: 


7T 

a^  =  1/7T  /  [6(t)-6(t+7r)  -  6(t-7i)]  cos  kt  dt  (111) 

— 7T 

With  the  help  of  Equation  21. 9- 10a  from  [31],  Equation  111 
may  be  integrated  to  give: 

afc  =  0  ,  k  -  even 

ak  =  2/t r,  k  -  odd  (112) 

As  was  done  for  PQ,  the  series  in  Equation  109  can 
be  rewritten  in  terms  of  the  delta  function  and  an  arbi¬ 
trary  function  H(<pi+cp)  as  follows: 

i 

=  cos  (p[hQ  -  (4/7r-2h0)  cos  2  C qpi+cp)  + 

(4/7r-2hQ)  (tt/2)  6 (cp^+q?)  +  H(cpi+cp)  ]  (113) 

where  H(qpi+qp)  is  subject  to  the  conditions: 

H(cp)  =  H(-tp) 

H(<p)  =  H(7T-<p) 
tt/2 

J  cos^  <p  H(tjH-cp^)  dqp  =  0  (114) 

-7T/2 
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Equation  113  satisfies  Equation  66-c  for  arbitrary  value  of 
tig  and  H  subject  to  the  conditions  set  forth  in  Equation  114 
and  such  that  P^  is  non-negative.  It  will  be  shown  in 
Section  2.4.3  that  tig  and  H  have  additional  restrictions 
on  them. 

2.4  INVARIANCE  OF  KERNEL  UNDER  ROTATION 
OF  COORDINATE  SYSTEM 

Any  mathematical  description  of  the  scattering 
kernel  must  be  of  the  form  such  that  P(v>t  ,~v)  is  invariant 
with  respect  to  a  rotation  of  the  coordinate  system  about 
the  surface  normal.  This  condition,  like  the  non-nega¬ 
tivity  of  the  scattering  kernel,  is  one  that  is  difficult 
to  apply  a  priori  but  is  simple  to  verify  after  a  solution 
is  obtained.  The  scattering  kernel  as  given  by  Equations 
13,  51  and  64  must  therefore  be  invariant  as  the  coordinate 
systems  are  rotated  (e.g.,  as  ^  is  varied).  Referring  to 
Figure  1,  Appendix  A,  it  can  be  seen  that  as  the  x,  z  axes 
are  rotated  about  the  y  axis  that  v  ' ,  v  ' ,  v  ,  v  ,  y .  , 

6,  <pi,  and  cp  vary,  while  vy ' ,  vy,  vrf,  vr>  v  ' ,  and  v 
are  invariant.  In  order  to  prove  that  a  scattering  kernel 
satisfies  this  condition  it  is  necessary  to  show  that  it  is 
invariant  with  ■•}/  for  all  possible  incident  and  reflected 
velocities.  On  the  other  hand  a  scattering  kernel  can  be 
disqualified  by  demonstrating  one  example  that  fails  to  hold. 
Each  of  the  three  kernels  in  each  coordinate  system  will  be 
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scrutinized  in  order  to  eliminate  those  that  are  unsatis¬ 
factory  and  thereby  simplify  the  equations. 


2.4.1  Check  of  Equations  in  Rectangular  Coordinates 

The  complete  scattering  kernel  is  given  by  Equation 
13  with  the  coordinate  kernels  given  by  Equation  27.  Remem¬ 
ber  that  the  form  for  the  z  coordinate  is  identical  with 
that  for  the  x  coordinate.  P(v',v)  must  be  independent  of 
the  orientation  of  the  x,  z  axes  for  arbitrary  values  of  An 
and,  in  particular,  for  An=l  with  all  other  A's  equal  to 
zero. 

For  n=0,  in  the  x  and  z  coordinates,  using  Equa¬ 
tion  38  for  P  n  and  P  n,  Equation  13  becomes: 

X  y  U  Z  ,  U 


Py(7r/2RTw)AxAz  exp 


{-[ 


,  1/2 

(ax2-i>  sx'-axsx 


[ 


(Az2-1) 


1/2 


S  '-A  S 
z  z  z 


n 


(115) 


The  expression  for  Py  was  not  required  since  it  is  a  con¬ 
stant  with  the  rotation  of  the  coordinate  system  since  vy 
and  v  *  do  not  change.  Equation  115  can  be  written  in  terms 
of  the  single  independent  variable  f  by  the  use  of  Table  I, 
page  18,  which  gives: 


P  =  7t/2RT  P  „A  A  exp 
w  y,0  x  z  * 


{-['Ax2- 


1/2  2 

1)  S  '  cos(Q+tf)  -  A  S  cos’!'] 
r  x  i* 


2  1/2  2\ 
[(Az  -1)  Sr’sin(n+T^)  -  AzSrsinf  ]  \ 


(116) 
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where  D,  = 

Notice  that  Q  is  independent  of  the  rotation  of  the  axis 
since  both  and  ^  vary  in  the  same  manner. 

The  scattering  kernel  given  by  Equation  116  is  in¬ 
variant  with  $  only  if  the  coefficients  of  terms  containing 
i/  vanish.  Expanding  the  argument  of  the  exponential  func¬ 
tion  in  Equation  116,  collecting  terms  with  like  variations 
of  f,  and  equating  the  coefficients  to  zero,  it  can  be 
shown  that  Equation  116  is  invariant  with  f  if  and  only  if 
A  =  A  .  For  this  condition  Equation  116  becomes: 

X  Z 

P  =  <tt/2RTw)  PyAx2  exp[-(Ax2-l)  Sr'2  + 

2Sr'SrAx(Ax2-l)1/2  cos  Cl  -Ax2Sr2]  (117) 

Thus  Equation  116  is  a  valid  scattering  kernel  only  for  an 
isotropic  surface. 

For  n  greater  than  zero  it  is  easy  to  show  that  the 
scattering  kernel  from  Equation  13,  using  the  numerical 
results  from  Figures  4  and  5,  Appendix  A,  is  not  invariant 
with  rotation.  Looking  at  the  results  for  n  =  2  and  4  it 
can  be  seen  that  for  each  value  of  a  ,  P  (and  therefore  P  ) 

X  X  z 

is  equal  to  zero  at  one  or  more  values  of  Sx  for  each  value 
of  S'.  A  consequence  of  the  assumed  product  solution  of 

A 

P  is  that  the  zeroes  form  an  intersecting  grid  of  infinite 

straight  lines  in  the  v  ,  v  (or  S  ,  S  )  plane.  It  is  easy 

x  z  x  y 

to  see  that  as  the  coordinate  system  is  rotated  the  zeroes 
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of  P  turn  also  and  remain  orthogonal  to  the  axis,  therefore 

P  is  not  invariant  with  rotation  for  at  least  n  =  2  and  4. 

It  is  assumed  that  P  for  n  greater  than  4  is  also  an  un- 

x,  n 

suitable  kernel. 

It  can  be  concluded  that  AQ  must  be  equal  to  unity 

with  the  other  A  *s  equal  to  zero,  while  B„  is  still  arbi- 

n  n 

trary  subject  to  the  previously  stated  conditions.  The 

final  results  for  P  and  P  are  therefore  given  by  Equation 

X  z 

38  with  A  =  A  . 

x  z 

2.4.2  Check  of  Equations  in  Cylindrical  Coordinates 

Equation  51  gives  the  scattering  kernel  in  cylin¬ 
drical  coordinates  with  the  individual  kernels  given  by 
Equations  27-b,  59  and  63.  It  is  obvious  that  Pr,  like  Py, 

is  invariant  with  rotation  of  the  x,  z  axes  since  v„  and 

r 

vr’  do  not  change.  From  Equation  63  is  equal  to  a  con¬ 
stant  and  is  consequently  independent  of  the  orientation. 
Therefore,  Equation  51  as  derived  satisfies  the  invariance 
with  rotation  condition. 

2.4.3  Check  of  Equations  in  Spherical  Coordinates 

The  scattering  kernel  in  spherical  coordinate  is 
given  by  Equation  64.  Because  the  speeds  v'  and  v  are  not 
affected  by  the  coordinate  system,  P^  is  invariant.  Using 
Equations  99  and  113  in  Equation  64  yields: 

P  =  Pp  cos  6  cos  cp[(l-e1)/2  +  e^  cos  (£^+0)  -e^  7r/2  6(0^+0)  + 
G(@i+e)][h0  -  (4/7r-2h0)  cos  2(cpi+cp)  + 
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(4/7T-2h0)  (7T/2)  6^+qp)  +  H(cpi+cp)  ]  (118) 

From  Figure  1,  Appendix  A,  and  Table  I,  page  18,  it  can  be 
seen  that: 


cos  0  cos 


(119) 


which  is  invariant  with  changes  in  iK  Therefore,  the 
product  of  the  two  bracket  terms  in  Equation  118  must  be 
invariant  with  axis  rotation. 

It  will  now  be  shown  that  e^^  and  4/7r-2hQ  .must  be 
zero  in  order  to  maintain  the  invariance  of  P.  The  term 
6(0i+0)  is  infinite  along  the  arc  0  =  -0^  and  6(qp^+qp)  is 
infinite  along  the  arc  qp  =  -qp^.  Since  neither  of  the 
bracketed  terms  is  identically  equal  to  zero  over  any  range 
of  0  or  qp,  P  is  infinite  (except  possibly  for  isolated 
points)  along  0  =  -0i  and  qp  =  -cp^.  (Note  that  being  in¬ 
finite  does  not  disqualify  a  term.)  However,  as  the  x,  z 
axes  are  rotated,  the  loci  of  constant  0  and  qp  rotate  about 
the  point  where  0  =  -0^  and  qp  =  -qp^.  This  means  that  the 
curves  where  P  is  infinite  moves  in  phase  space  and  there¬ 
fore  P  cannot  be  invariant  unless  the  coefficients  of  the 
delta  function  are  zero.  This  means  that: 


el 


=  0 


(120) 


and , 


h 


0 


2 /it 


(121) 
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In  a  similar  manner  it  will  be  shown  that  G  and  H 
in  Equation  118  must  be  identically  zero.  From  Equation 
100  it  can  be  seen  that: 

G(tt/2)  =  -G(tt/2)  =  0  (122) 

Since  e^  is  zero  from  Equation  120  and  from  Equation  122 
0(0+0^  is  zero  for  0+0^^  =  7t/2  the  first  bracket  term  in 
Equation  118  becomes  1/2  for  0  =  it/2-9^.  Since  it  has  been 
shown  that  the  product  of  the  two  bracketed  terms  must  be 
a  constant,  the  second  term  must  be  a  constant  for  arbi¬ 
trary  value  of  <p  and  The  only  solution  for  G  and  H  is 

therefore: 


G  -  0 
H  3  0 


(123) 


Equations  99  and  113  become,  respectively: 

pe  "  ££f-2  <124) 

,  2  s°5-2.  (125) 

<P  -  7T 

The  scattering  kernel  in  spherical  coordinates  that 
satisfies  all  necessary  conditions  is  given  by  Equations  64, 
74,  124  and  125.  For  the  case  of  Dn  =  1,  D^=0  i/n  and 
Ap— »*oo  it  can  be  shown  that  the  scattering  kernel  in  spher¬ 
ical  coordinates  approaches  that  for  diffuse  elastic  scat¬ 
tering  as  given  by  Equation  9. 


58 


A  E  DC -TR -7  2-1 58 


III.  APPLICATION  OF  THEORETICAL  RESULTS 

In  Section  I  the  most  general  description  of  gas- 
surface  interaction,  that  of  the  scattering  kernel,  was 
introduced  and  in  Section  II  several  scattering  kernels  were 
developed  that  satisfy  all  required  conditions  including  that 
of  reciprocity.  However,  the  scattering  kernel  itself  is 
far  too  detailed  to  be  of  much  practical  use  since  it  gives 
the  probability  density  as  a  function  of  six  independent 
velocity  components  as  well  as  characteristics  of  the  sur¬ 
face.  Usually  it  is  some  moment  of  the  reflected  molecular 
distribution  such  as  the  molecular  flux  distribution,  mean 
velocities,  or  total  energy  that  is  desired.  In  this  section 
several  of  these  important  moments  will  be  obtained,  some  in 
closed  mathematical  form  and  others  from  numerical  integra¬ 
tion. 


3.1  REFLECTED  GLOBAL  PROPERTIES 

Mean  reflected  molecular  properties  are  determined 
by  integrating  the  desired  properties  over  some  group  of 
reflected  molecules.  Global  properties  are  determined  by 
integrating  over  all  reflected  molecules.  Total  energy  and 
momentum  exchange  between  a  gas  and  a  surface  are  important 
quantities  which  are  determined  in  the  remainder  of  this 
section . 
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3.1.1  Mean  Reflected  Velocity  Components 

In  order  to  determine  the  momentum  exchange  between 
a  gas  and  a  surface  it  is  necessary  to  know  the  mean  veloc¬ 
ities  of  reflected  molecules.  The  mean  velocity  component 
in  a  given  direction,  say  in  the  direction  of  a  unit  vector 
k~  of  reemitted  molecules  is  given  by: 

v-k  =  V  =  /  TT*\Tp  d3v  (126) 

v*n>0 

3.1. 1.1  The  Determination  of  Sx 

In  rectangular  coordinates,  using  Table  I,  page  18, 
and  Equation  126,  the  mean  velocity  component  in  the  direc¬ 
tion  of  the  x  coordinate,  v  ,  is: 

X 

00  03  00 

vx  -  /  /  /  v  P  dv  dv  dv  (127) 

-oo  0  3 

Since  it  has  been  assumed  that  P  is  given  as  a  product  solu¬ 
tion  in  Equation  13,  Equation  127  can  be  separated  into  the 
product  of  three  integrals  as  follows: 


00 


oc 


v  P  dv 

XX  X 


/ 

o 


P  dv  / 
v  V  J 
J  J  -00 


Vvz 


(128) 


It  can  be  seen  that  the  last  two  of  these  integrals  are 
unity  from  Equations  14-b  and  14-c.  Using  this  fact  and 
the  result  of  Section  2.4.1  that  P  is  given  by  Equation  38, 
Equation  128  becomes  in  terms  of  dimensionless  velocities: 
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oo 


Sx  = 


-7=.  / 

V7T  —00 


Sv  exp  [- 


(A  2-l)1/2 
x 


S  ’-A  S 

X  XX 


]  dS 


x 


(129) 


Integrating  yields: 


or  in  terms  of  ax : 


(130) 


(131) 


Thus,  the  mean  reflected  x  component  of  velocity  is  propor¬ 
tional  to  the  x  component  of  the  incident  velocity. 

In  both  cylindrical  and  spherical  coordinates,  S 

A 

is  zero  from  symmetry. 


3. 1.1. 2  The  Determination  of  Sy 

Again  using  Table  I,  page  18,  and  Equation  126,  the 
mean  velocity  component  in  the  y  direction,  vy,  is: 

00  00  on 

V  =  /  f  /  v  Pdv  dv  dv  (132) 

y  -oo  o  -oo  y  x  y 


for  the  rectangular  coordinate  solution.  As  was  the  case 
for  v  ,  Equation  132  reduces  to: 


=  / 


v  P  dv 

y  y  y 


(133) 


Where  Py  is  given  by  Equation  27-b.  Equation  133  can  be 
evaluated  for  each  P  .  For  n=0,  using  Equation  47, 

y  * n 
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Equation  133  becomes  in  terms  of  dimensionless  velocities: 

Sy  ‘  ay2  l  sy2  Vsy]  a*p[-<Ay2-i)Sy'2  - 

Ay2  sy2]  dSy  (134) 

Using  the  results  of  Appendix  B,  Equation  134  becomes: 

Sv  =  exp[-(A  2-l)S'2]  M[3/2,  1,  (A  2-1)S  ’]  (135) 

y  2  a  y  y  y  y 

where  M  is  the  confluent  hypergeometric  function  of  Kummer. 

The  result  of  Equation  135  is  shown  in  Figure  14, 
Appendix  A,  with  S  given  as  a  function  of  S  *  for  various 

y  y 

values  of  ay. 

For  n  greater  than  zero  it  was  necessary  to  evaluate 
Equation  133  numerically.  Using  the  numerical  results  for 
Py  n  obtained  in  Section  2.1.2,  Equation  133  was  integrated. 
The  results  for  n=l  and  2  are  shown  in  Figures  15  and  16, 
Appendix  A,  respectively,  for  given  values  of  ay. 

The  result  of  in  cylindrical  coordinates  is 
exactly  the  same  as  that  given  above  since  it  was  shown  in 
Section  2.2  that  P  is  common  to  both. 

y 

The  solution  for  P  in  spherical  coordinates  from 
Equation  64  yields  for  the  mean  normal  velocity: 

00  tt/2  tt/2 

W  ;  /  v^  cos2  <p  cos  Q  PpP^P^dOdcpdv  (136) 

0  -7T/2  -7T/2 
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3 

Use  of  Table  I,  page  18,  was  required  to  express  vy  and  d  v 
in  spherical  coordinates. 

Separating  variables  and  substituting  Equations  124 

and  125  for  P-  and  P  ,  respectively,  Equation  136  becomes: 

y  (p 


v 


y 


o 


dv 


K  2, 

7T/,2  ,  2 

f  003  8  d0 

r  2  cos  cp 

J  2 

J  * 

-7T/2 

-tt/2 

dtp 


(137) 


Where  Pp  is  given  by  Equation  74  and  vy  must  be  evaluated 
for  n=0, 1,2,  etc. 

For  n=0  (D  =1),  using  Equation  83  for  P^  n  and 
n  p ,  u 

performing  the  two  angular  integrations,  Equation  137  gives 
for  the  dimensionless  mean  reflected  normal  component: 


(-S’) 


/  S3  exp[-(A  2-l)S'2  - 
0  ’  H 


A 


P 


2 


IX[-2S’S 


(138) 


Using  the  result  of  Appendix  B,  Equation  138  becomes, 


Sy  =  2A^  exP[-CAp2-1)S'2]M[5/2,2,S’2(Ap2-l)]  (139) 

The  result  of  Equation  139  is  shown  in  Figure  17, 

Appendix  A,  with  plotted  as  a  function  of  S'  for  various 

values  of  a  .  The  numerical  results  for  n=l  and  2  are 
P 

shown  in  Figures  18  and  19,  Appendix  A,  respectively. 
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3.1.2  Mean  Energy  of  Reflected  Molecules 

The  exchange  of  energy  between  a  gas  and  a  surface 
is  perhaps  the  most  important,  or  at  least  the  most  investi¬ 
gated,  gas-surface  phenomenon.  The  kinetic  energy,  e,  of  a 
reemitted  molecule  is  given  by: 


e 


(140) 


The  mean  reflected  energy  is: 


e 


„/ 

v  n<0 


e  P  d^v  = 


®  v2 
2 


(141) 


The  evaluation  of  e  thus  requires  the  determination  of  the 
mean  square  of  the  reflected  velocities. 


3.1. 2.1  The  Determination  of  S2  in  Rectangular  Coordinates 

In  rectangular  coordinates  the  square  of  the  velocity 
is  given  by: 


2 

v 


2  2 
V  +  V 

X  y 


(142) 


Substituting  this  expression  into  Equations  131  and  132  and 
using  Equation  5,  the  mean  square  of  the  velocity  of  re¬ 
flected  molecules  becomes: 


—  00 


00 


/ 

0 


/ 


—  00 


P  P  P  dv  dv  dv 
x  y  z  x  y  z 


(143) 


Separating  variables  and  utilizing  the  normalization 
conditions  (Equation  6)  Equation  134  becomes: 
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00  O  00  O  00  o 

v'2  -  f  v  dv  +  |  v  dv  +  f  v  1  dv 
xxx  iyyy  J  zzz 


0 


(144) 


which  is  equivalent  to: 


~7  — 2  — 2  — 9 

v  =  v  +  v  +  v  (145) 

x  y  z 

In  terms  of  dimensionless  velocities  the  terms  on 
the  right  hand  of  Equation  145  become: 

— 2  00  2 

S  -  /  S  2RT  P  dS  ( 146-a) 

x;x  wxx 

—  co 


2RT  P  dS 
w  y  y 


(146-b) 


— 2  2 
where  the  form  of  S  is  identical  to  that  of  S 

For  n=0,  Equations  29  and  38  may  be  utilized  for 
j2RTy  Px  and  ^2RT^  P^ ,  respectively,  in  Equation  146  and 
integrated  to  yield: 


(147-a) 


(147-b) 


Use  was  made  of  Appendix  B  for  obtaining  Equation  147-b. 
Figure  20,  Appendix  A.  presents  Equation  147-a  with 
plotted  as  a  function  of  S  '  for  various  values  of  a  .  Fig- 
ure  21,  Appendix  A,  similarly  shows  the  results  of  Equation 
147-b. 
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For  n  greater  than  unity  it  has  been  shown  (Section 
2.4.1)  that  An  =  0  and  therefore  the  complete  scattering  ker¬ 
nel  in  the  x  and  z  coordinates  is  given  by  Equation  38.  In 
the  y  direction,  the  numerical  results  for  n  greater  than 
unity  as  shown  in  Figures  6  and  7,  Appendix  A,  was  substi¬ 
tuted  in  Equation  146-b  and  numerically  integrated.  These 
results  for  n=l  and  2  are  shown  in  Figures  22  and  23,  Appendix 
A. 


3. 1.2.2  The  Determination  of  S2  in  Cylindrical  Coordinates 

In  cylindrical  coordinates  the  square  of  the  molecu¬ 
lar  velocity  is  given  by: 

22  2 

v  =  vr  +  vy  (148) 

Substituting  Equation  148  into  Equations  140  and  141 
and  using  Equation  51  yields: 


00 


/ 

0 


00 


/ 

0 


2  2 

(v  +v  )  P  P  P.v  dv  dv  df 
r  y  y  r  if  v  y  r 


(149) 


After  separation  of  variables  and  the  utilization 
of  the  normalization  conditions  as  expressed  in  Equation  53, 
Equation  149  becomes: 


d  v 


r 


00 


+  / 
0 


which  is  the  same  as: 


(150) 
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2  2,2  n,n 

v  =  vr  +  vy  (151) 

2 

The  results  for  vy  are  identical  to  those  in  rectangular 
coordinates  and  already  presented  in  Figures  21,  22  and  23, 
Appendix  A,  as 

In  terms  of  dimensionless  velocities  the  first  term 
on  the  right  hand  side  of  Equation  151,  after  using  Equa¬ 
tion  62  and  integrating,  becomes  for  n=0: 


+  S 


r 


(152) 


which  is  of  the  same  form  as  Equation  147-b.  In  fact,  from 
an  inspection  of  Equations  61,  146-b,  and  the  expression  for 
vr2,  it  can  be  seen  that  the  results  for  vy  and  vr  are 
equivalent  for  all  values  of  n.  Thus  Figures  21,  22  and 
23,  Appendix  A,  also  give  the  results  for 


3.1. 2.3  The  Determination  of  S2  in  Spherical  Coordinates 

From  Table  I,  page  18,  it  can  be  seen  that  the 
square  of  the  molecular  velocity  in  spherical  coordinates 
is  simply: 


v  =  vp*  (153) 

Equation  141,  after  integrating  over  9  and  cp  becomes: 


P  dv 


(154) 
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For  n=0,  using  Equation  83  and  Appendix  B,  Equation  154 
becomes  in  terms  of  dimensionless  velocity: 


(155) 


S2  =  2A  “2  +  S’ 2(1-A  ~2) 

P  P 


Equation  155  is  plotted  in  Figure  24,  Appendix  A,  as  J  S 


as  a  function  of  S'  for  various  values  of  a  .  For  n  greater 

P 

than  1  it  was  necessary  to  numerically  integrate  Equation 
154  using  the  numerical  results  for  shown  in  Figures 
12  and  13,  Appendix  A.  The  results  are  presented,  using 
the  same  parameters  as  for  n=0,  in  Figures  25  and  26,  Ap¬ 
pendix  A,  for  n=l  and  2,  respectively. 


3.2  REFLECTED  DIRECTIONAL  PROPERTIES 


In  Section  3.1  the  reflected  global  properties  of 
velocity  components  and  energy  were  determined  based 
upon  the  theoretical  results  of  Section  II.  We  will  now 
turn  our  attention  to  the  determination  of  local  or 
directional  properties. 

Consider  the  detector  shown  in  Figure  27,  Appendix 
A.  The  type  of  detector  is  immaterial,  it  could  measure 
molecular  flux,  density,  velocity,  etc.  The  detector, 
located  a  distance  L  from  the  point  of  impact,  Q,  of  the 
molecular  beam,  has  a  sensing  area,  dA,  which  subtends 
a  solid  angle  of  do>  when  viewed  from  A.  The  solid  angle 
is  given  by: 
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du  =  ^  (156) 

L2 

Because  all  reflected  molecules  are  assumed  to 
emanate  from  the  origin,  all  molecules  within  dau  at  9  and 
cp  in  physical  space  are  also  within  the  same  solid  angle -in 
velocity  space.  Thus,  the  molecular  properties  at  the 
detector  are  determined  by  integrating  the  desired  pro¬ 
perties  over  all  molecules  within  dco. 

Previous  results  in  this  paper  have  been  presented 
(e.g.,  Figures  2-26,  Appendix  A,  inclusive)  for  each  of 
three  coordinate  systems  for  various  values  of  n.  Inas¬ 
much  as  the  scattering  kernel  with  n=0  in  rectangular 
coordinates  embodies  most  of  the  experimentally  observed 
gas-surface  phenomena,  the  results  which  follow  will  be 
presented  in  detail  only  for  this  case,  since  more  in¬ 
dependent  variables  are  involved  because  the  coordinates 
directions  are  not  separable  as  was  true  for  the  global 
properties. 

Important  directional  properties  which  will  be 
investigated  are  molecular  flux,  density,  speed  and 
energy. 

3.2.1  Molecular  Flux  Distribution 

The  easiest  reflected  quantity  to  measure  is 
molecular  flux.  Referring  to  Figure  27,  Appendix  A,  the 
number  of  reflected  molecules  with  speed  between  v  and 
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v  +  dv  entering  the  detector  during  time  dt  will  now  be 
determined.  The  number  of  molecules  reflected  during 
time  dt  with  speeds  between  v  and  v  +  dv  within  doi  is: 

NP  v2  doi  dv  dt  (157) 

2 

Notice  that  v  dec  dv  is  simply  the  incremental  volume  in 

velocity  space.  These  molecules  are  spread  out  over  a 

length  vdt  within  da>  in  physical  space.  The  number  per 

unit  volume  in  physical  space  is  that  given  by  Equation 

2 

157  divided  by  the  increments  volume,  L  dcovdt,  which  yields: 


dp 


NPvd  v 


(158) 


During  the  time  dt,  molecules  with  velocity  v  contained  in 
the  volume  dA  vdt,  will  enter  the  detector.  This  volume 
multiplied  by  the  incremental  density  given  by  Equation  158 
and  dividing  by  dt  in  order  to  obtain  an  intensity  yields: 


dl 


NP  v  dvdA 


(159) 


dl  is  therefore  the  total  number  of  molecules  with 
speeds  between  v  and  v  +  dv  that  enter  the  detector  per 
unit  time.  Equation  159  may  be  integrated  over  all  re¬ 
flected  speeds  to  obtain  the  total  number  of  molecules 
scattered  into  the  detector  per  unit  time.  Thus: 
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Combining  Equations  160  and  161  yields  in  terms  of  dimen¬ 
sionless  speed: 


I  =  (2RT  )3/2  /  PS2dS  (162) 

W  0 

From  a  consideration  of  the  normalization  condition 
given  by  Equation  2  it  can  be  shown  that  I  is  the  proba¬ 
bility  density  that  a  reflected  molecule  will  be  in  a  given 
direction  without  regard  to  speed.  Thus  I  den  is  the  pro¬ 
bability  that  a  molecule  will  be  reflected  within  den,  and 
therefore  enter  the  detector. 

Using  Equation  13  for  the  scattering  kernel  with 
the  individual  kernels  given  by  Equations  38  and  47, 

Equation  162  becomes,  in  terms  of  dimensionless  velocities: 

9  9 

2A  ^  00  o  T  /o  o 

I  =  — - — —  /  expf-[(A  -1)L/  S  '-AS  cosqpsine] 

7T  q  v.  X  XX 

Ax2S2  sin2  cp  -(Ay2-l)Sy2-Ay2S2cos2  q>  cos2  0)  IQ[2Ay(Ay2- 

1)1//2  Sy’S  cos  cp  cos  0]  S3  cos  <p  cos  0  dS  (163) 
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Oj 

For  a  given  A  and  A  ,  I  can  be  evaluated  numerical- 
x  y 

ly  from  Equation  163  as  a  function  of  cp  and  0  for  various  com¬ 
binations  of  S  *  and  S  1 .  (Integration  in  closed  form  was 

x  y 

not  achieved.)  The  numerical  results  of  Equation  163  for 

a  =  0.5  and  a  =  0.8  are  shown  in  Figure  28,  Appendix  A, 
x  y 

for  0^  =  -20°,  -50°  and  -80°.  Similar  results  for  ax  =  0.8 
and  ay  =  0.5  are  shown  in  Figure  29,  Appendix  A.  The  numer¬ 
ical  results  are  presented  as  a  function  of  0  for  qp  =0,  and  as  a 

function  of  ©  for  0=0  where  0  is  the  value  of  0  for  which 

m  m 

I  is  a  maximum.  The  results  in  Figures  28  and  29  are  typical 

for  a  <  a  and  a  >  a  ,  respectively, 
x  y  y  x 

For  the  special  case  of  A  =A  =A(a  =  a  )  Equa- 

x  y  x  y 

tion  163  may  be  simplified  somewhat.  For  this  condition 
Equation  163  becomes: 


I  =  —  f  exp[-(A2-l)S’2  +  2A./a2-1  S'  sin  0.  S  cos  qp  - 
77  0  V 

^2g2]  Iq[2A  (A2-1)1'^2S,S  cos  0^  cos  cp  cos  0]  • 

(AS)  2  cos  qp  cos  0  d(AS)  (164) 


Equation  164  may  be  integrated  numerically  to  ob- 


tain  I  as  a  function  of  0r>  0^  and 


The  results 


of  Equation  164  for  cp  =  0  and  0=0  are  shown  as  a  function 

_____  HI 

of  0,  and  ^A2-l  S’  for  0£  of  -20°,  -50°  and  -80°  in  Figure 
30,  Appendix  A. 
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3.2.2  Molecular  Density  Distribution 

The  total  density  of  reflected  molecules,  p,  is  ob 
tained  by  integrating  Equation  158  over  all  speeds.  Thus: 


N 


00 


P  =  “9  /  P  v  d  v 
0 


(165) 


A  dimensionless  density,  p,  is  defined  as  follows: 
.  2 

Jf  =  - 

N  2RT 

w 


(166) 


Combining  Equations  165  and  166  gives  for  p  in  terms  of 
dimensionless  speed: 


p  =  (2RT  )3/2  /  P  S  d  S  (167) 

W  0 

The  numerical  results  for  p  are  very  similar  to 
those  for  I*  and  therefore  will  not  be  presented  at  this 
time.  However,  Equation  167  will  be  utilized  later  when 
comparison  of  theory  to  experimental  data  is  made. 

3.2.3  Mean  Reflected  Speed 

Mean  reflected  speed  is  defined  as  the  average  speed 
of  molecules  entering  the  detector;  i.e.,  the  average  is 
over  the  flux  of  molecules  and  not  over  density  as  is  usually 
the  case  in  the  kinetic  theory  sense.  Using  Equation  159 
for  the  incremental  flux  of  molecules,  the  mean  velocity 
becomes: 
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“  3 

v  =  (I  )  /  Pv  dv  (168) 

0 

or  in  terms  of  dimensionless  velocity: 

S  =  (l“1)(2RT  )3/2  /  P  S3  d  S  (169) 

W  0 

Using  Equations  13,  38  and  47,  Equation  169  may 
be  written  in  a  form  similar  to  Equation  163  with  I  evalu¬ 
ated  from  Equation  163.  Numerical  results  for  S  for  ax  = 

0.5  and  ay  =  0.8  are  shown  in  Figure  31,  Appendix  A,  for 

0.  =  -20°,  -50°  and  -80°  for  cp  =  0  and  8=6.  Similar 
l  m 

results  for  ax  =  0.8  and  ay  =  0.5  are  shown  in  Figure  32, 

Appendix  A.  The  results  for  S  presented  in  Figures  31  and 

32  are  typical  for  a  <  a  and  a  >  a  ,  respectively. 

x  y  y  x 

For  the  special  case  of  A  =  A  =  A  the  number  of 

x  y 

independent  variables  in  Equation  169  may  be  reduced  to 

0,9,0^  and  \/a2-1  S'  with  AS  as  the  dependent  variable. 

Figure  33,  Appendix  A,  gives  the  numerical  results  for 

AS  with  0.  =  -20°,  -50°  and  -80°  for  qp  =  0  and  0=0. 

l  m 

3.2.4  Mean  Reflected  Energy 

Mean  reflected  energy  like  mean  speed  is  determined 
by  averaging  over  the  flux  of  molecules.  Using  Equation 
140  the  mean  energy  of  molecules  reflected  in  a  given 
direction  is : 

e  =  (I-1)  /  P  v2  dv  (170) 
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which  is  equivalent  to: 


e 


m  7 


or  in  terms  of  dimensionless  velocities: 


(171) 


e 


=  ~  2RT  ? 
2  w 


(172) 


where: 


?  =  (I_1)  /  (2RT  )3/2  P  S4  d  S 
0  W 


(173) 


The  numerical  results  for  y?  are  very  similar  to 
those  for  S  and  will  therefore  not  be  given  as  such.  However, 
the  directional  energy  accommodation  coefficient  a  is  often 
used  to  describe  energy  exchange  with  the  surface.  It  is 
defined  in  terms  of  molecular  energy  as  follows: 


e  -  e. 

a  -  - ±- 

e  e„.  -  e  • 


(174) 


where  e^  is  the  molecular  incident  energy  and  ew  is  the 
reflected  energy  if  at  equilibrium  with  the  surface  at  Tw- 
Equation  174  may  be  rewritten  in  terms  of  dimensionless 
speeds  as  follows: 


(175) 
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Use  was  made  of  the  fact  that  e  =  2mRT  . 

w  w 

The  local  value  of  a  is  not  bounded  by  zero  and 

e 

unity  as  one  normally  expects.  From  the  second  law  of 
thermodynamics  the  total  energy  accommodation  must  be  within 
these  bounds.  However,  there  are  no  such  restrictions  upon 
the  local  value.  In  fact,  from  Equation  175  it  can  be  seen 
that  for  S.  =  \[2,  a  diverges  unless  =  S.^,  which  occurs 

1  6  X 

only  in  isolated  directions.  This  uncertainty  raises  serious 
questions  concerning  the  usefulness  of  a0,  however,  it  is 
widely  used  in  the  literature  and  it  is  introduced  here  for 
that  reason. 

3.2.5  Direction  of  Maximum  Molecular  Flux 

The  direction  of  maximum  molecular  flux  is  one  of 
the  most  widely  used  parameters  in  the  study  of  gas-surface 
interactions.  It  seems  that  the  test  of  any  gas-surface 
model  or  theory  is  to  compare  the  predicted  maximum  in  re¬ 
flected  flux  to  experimental  data  or  to  other  theories 
(e.g.,  see  the  review  paper  by  Trilling  [32]). 

These  maxima  in  flux  for  the  present  theory  can  be 

determined  for  typical  values  of  a  ,  a  and  0.  from  the  flux 

x  y  i 

distributions  as  presented  in  Figures  28,  29  and  30,  Ap¬ 
pendix  A.  Numerically  differentiating  these  results,  the 
value  of  9  for  which  I  is  a  maximum  can  be  determined,  this 

9  value  will  be  denoted  by  9  .  Value  of  0  obtained  for 

m  m 

the  values  of  a  ,  a  and  9.  given  previously  are  presented 
x  y  i 

in  Figures  34,  35  and  36,  Appendix  A. 
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The  limit  of  9  for  large  values  of  S'  can  be 
m 

determined  in  closed  form.  Recalling  that  in  Sections 

2.1.1  and  2.1.2  it  was  shown  that  as  A  and  A  increase 

x  y 

without  limit  specular  reflection  is  approached.  In  a 

similar  manner  as  S  '  and  S  1  increase  without  limit, 

x  y 

the  scattering  kernels  given  by  Equations  29  and  38  ap¬ 
proach  the  following: 

px— “  57  6 1  AxWj  •  V— ”  <176> 

^  py—  5^-  6  [f7^+  \vv] '  v— ' •  (177) 


Substituting  into  Equations  128  and  133,  S  and  become: 

x  y 


(178) 


Sy  '  ^7Sy 


S  ' 

y 


(179) 


Since  for  S '  —  oo  both  P  and  P  approach  a  delta 

x  y 

function  the  flux  of  molecules  also  approaches  a  delta 

function,  i.e.,  a  spike  directed  along  the  direction 

determined  by  S~~  and  S- .  Thus,  9  is  given  by: 

x  y  m 


8m  =  Tan"1  (Sx/Sy),  S’- 


(180) 


or  using  Equations  178  and  179  and  writing  in  terms  of  6K 
yields; 
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9  =  Tan”1  [-(Jl-a/Jl-a)  Tan  0.]  ,  S’ — ► »  (181) 

ni  x  y  l 

Figure  37,  Appendix  A,  presents  the  results  of 

Equation  181.  Notice  that  only  for  the  case  of  a  =  a  is 

x  y 

the  specular  limit  approached  as  S  r  00  and  depending  on 

a  and  a  the  limit  can  be  supra-specular  or  subspecular. 
x  y 

3.3  REFLECTED  VELOCITY  DISTRIBUTION  FUNCTION 

Throughout  Section  II  and  thus  far  in  Section  III 
we  have  considered  only  the  reflected  probability  of  a 
single  incident  molecule  or  the  mean  reflected  properties 
of  an  incident  monoenergetic  molecular  beam,  i.e.,  of 
uniform  speed  and  direction.  The  special  case  of  a  gas 
with  a  Maxwellian  distribution  incident  upon  a  surface 
will  now  be  considered  in  each  of  the  three  coordinate 
systems.  It  is  assumed  that  the  velocity  distributions  of 
incident  molecules  is  given  by: 

f"  =  (27rRT")"3/2  exp  (-v1 2/2RT”)  (182) 

The  velocity  distribution  function  f”  and  f+  were 
discussed  in  Section  1.2  and  the  relation  between  them  and 
the  scattering  kernel  was  given  in  Equation  3.  The  deter¬ 
minations  of  f+  will  be  made  only  for  those  scattering 
kernels  that  have  been  obtained  in  closed  form,  i.e.,  for 
n  =  0.  It  will  be  shown  that  f+  is  Maxwellian  in  form  in 
each  coordinate  system,  and  for  equality  of  accommodation 
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coefficients  thus  satisfies  the  so  called  "Knudsen's 
Assumption"  discussed  in  Section  1.1,  page  4. 


3.3.1  f*  in  Rectangular  Coordinates 

In  rectangular  coordinates,  using  Table  I,  page  18, 
and  Equations  13,  38,  47  and  182,  Equation  3  becomes: 


f 


+ 


where: 


f+  f+  f+ 
x  y  z 


(183) 


A  00 

f+  =  _ * 

X  TT  yj2RT~  —oo 


-  7  exp  [-(/a  2-l  S  ’-A  S  )2-T  S  ' 2  n~  ]  dS  ' 

—  j  t-LNVx  x  XX  w  x  J  x 


(184) 


,  ^A  0  /  o  p  p 

fy  =  £.1  _  f  Sy,I0(-2Ay^/Ay2-l  Sy'Sy)  exp[-Sy ,2(Ay  -1  + 


\/27rRT_ 


T  /T”)  -  A  2S  2]  dS  ' 

w  y  y  y 


(185) 


and  f+  is  of  exactly  the  same  form  as  Equation  184  with  S 

Z  X 

replaced  by  Sz- 

Integrating  Equations  184  and  185  (use  of  Appendix 
B  was  required  for  the  latter)  yields: 


_  -1/2 

,  A  (2ttRT  )  ,92  2  -  > 

ft. - ^ - t  7rT  exp{-SxZTwAx7[T  <AX  -1  +  T/T  >  0 


X  (A  2-l  +  T  /T")1/2 
'  x  w  7 


(186) 


f  = 


A  2(2ttRT")"1/2 

■  r 


(A  2-l  +  T  /T") 
v  y  w 


—  exp  f-S  2T  A  2/[t"(A  2-l  +■  T  /T”)  ]) 
i—  \  ^  \  y  w  y  1  y  w  V 


(187) 
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f2  =  fX  <SZ)  (188) 

Notice  that  the  reflected  velocity  distribution  is 

Maxwellian  in  form  with  possible  differences  in  temperature 

in  the  normal  and  tangential  directions.  Defining  the 

remitted  temperature  in  the  normal  and  tangential  directions 

as  T*  and  T+,  respectively,  Equations  186,  187  and  188 

y  * 

become: 


f+  =  (27rRTj)"1/2  exp(-vx2/2RT^)  (189) 

fy  =  (T"/Ty)1/2  (27rRTy)_1/2  exp(-vy2/2RTy)  (190) 

f*  =  (27rRT^)“1/2  exp(-vz2/2RTx)  (191) 

Solving  for  T*  and  T*  from  Equations  186,  187, 
x  y 

189  and  190  yields: 

Tx  =  T"(Ax2-1  +  VT")/Ax2  (192) 

Ty  =  T”  (Ay2-1  +  Tw/T_)/Ay2  (193) 

-2  -2 

Writing  Equations  192  and  193  in  terms  of  A  and  A  , 

x  y 

respectively,  yields: 

-2  Tx"T" 

*  -  — — r  (194) 

x  T  -T 


80 


A  E  DC -TR -72-1 58 


T+-t" 


T  -T 
W 


(195) 


Equations  194  and  195  are  the  normal  definitions 
of  the  thermal  accommodation  coefficients,  a  , .a  ,  general- 
ized  to  each  coordinate  direction.  Thus,  the  relation 
between  the  a's  and  A's  is  given  by: 


a  =  A 

x  x 


(196) 


a  =  A 
y  y 


(197) 


These  definitions  of  a  and  a  were  introduced  in  Section  II 

x  y 

without  proof. 

“  +  1  /2 

The  term  (T  /T  )  in  Equation  190  came  about  from 

y  y 

the  conservation  of  molecules  at  the  surface.  Thus,  the 
density  of  reflected  molecules  is  (T  /T  )  times  that  of 

y  y 

the  incident  density. 

For  the  special  case  of  A  =  A  (a  =a  =a)  Equation 

x  y  x  y 

183  becomes 


f+(v)  =  (T"/T+)  1/2  (27TRT+)"3/2  exp[-v2/2RT+]  (198) 


and  a  is  the  ordinary  thermal  accommodation  coefficient: 


T+-T~ 

T  -T” 
w 


(199) 
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3.3.2  t*  in  Cylindrical  Coordinates 

The  reflected  velocity  distribution  function  will  be 
evaluated  only  for  n  =  0.  In  cylindrical  coordinates,  using 
Table  I,  page  18,  and  Equations  51,  62,  47,  63  and  182, 
Equation  3  becomes: 

f+  =  f+  f+  (200) 

where  f*  is  given  by  Equation  185  and  after  integration  by 
Equation  190.  And: 

-  2  A^^RTT1  f  Sr'  I0<-2Aj/Ar2-l  Sr'Sr>- 

—  00 

exp[-Sr’2(Ar2-l  +  Tw/T")  -Ar2Sr2]  dS/  (201) 

Integrating,  using  Appendix  B  yields: 

f*  =  (27tRT^)_1  exp(-v2/2RT*)  (202) 


where: 


Tr  =  T  (Ar2-1  +  Tw/T‘)/Ar2 
_2 

or  solving  for  Ar  : 


A 


tJ-t 

T  -T 


(203) 


which  is  again  the  thermal  accommodation  coefficient,  ar« 
Thus: 
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ar  =  Ar-2  (204) 

For  the  special  case  of  Ar  =  Ay(ar=ay=a)  Equation 
200  becomes  identical  to  Equation  198. 

3.3.3  f+  in  Spherical  Coordinates 

In  spherical  coordinates  for  n  =  0  and  using  Table  I, 
page  18,  and  Equations  64,  83,  124  and  125,  Equation  3 
becomes,  after  integrating  out  the  angular  parts: 

f+  =  (2A  3/S)(A  2-l)"1/2(27rRT")"3/2  /  S'2  exp[-(A  2-l)S'2  - 

P  P  —00  P 

Ap2S2  -  TwS’2/T-)]  I1(-2S'SApyAp2-l)  dS'  (205) 

Using  Appendix  B  and  integrating  yields: 


f+(v)  =  (T"/T  1/2  (2ttRT+)"3/2  exp(-v2/2RT*) 

P  P  P 


where : 


T+  =  T  (A  2-l  +  T  /T")  A  2 
P  P  vr  p 


_2 

and  solving  for  Ap  gives: 


-2  =  T+-T 

3  T  -T~ 

w 


(206) 


(207) 


(208) 


which  is  again  the  thermal  accommodation  coefficient,  a  ■ 
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IV.  COMPARISON  OF  THEORY  WITH  EXPERIMENT 

The  purpose  of  this  section  is  not  to  define  general¬ 
ized  gas-surface  interacting  parameters,  but  to  show  that 
the  theoretical  results  developed  in  the  first  three  sections 
of  this  paper,  indeed  are  capable  of  describing  real  gas 
surface  effects.  Experimental  reflected  flux  distributions 
from  a  wide  range  of  gases  and  surfaces  indicate  that  certain 
scattering  characteristics  are  common  to  almost  all  condi¬ 
tions  . 


4.1  GENERAL  CHARACTERISTICS 

Stickney  [33,  page  186]  in  an  excellent  article  on 
scattering  from  solid  surfaces,  lists  eight  scattering 
characteristics  which,  in  terms  of  present  nomenclature,  are: 

1.  dem/dv’  >  0 

2.  d0  /dT  <  0  (This  characteristic  is  invalid  if 

m  w  — 

the  composition  or  structure  of  the  surface 

varies  with  T  .) 

w 

3.  aem/a(-ei)  >  o 

m  i  — 

4*  36_/ani  <  0  when  0  <  -0. 

m  mi 

d0  /dm  >  0  when  0  >  -0. 

m  mi 

(It  is  assumed  that  AH,  the  heat  of  adsorption, 
is  held  constant  in  these  partial  derivatives.) 

5.  The  dispersion  of  the  diffuse  scattering  patterns 
increases  with  AH  and/or  with  m. 
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6.  Supraspecular  patterns  (i.e.,  0m  >  -9^)  are 
most  likely  to  occur  when  both  of  the  following 
conditions  are  met:  (a)  S'  >  1  and  (b)  -0^  « 
90°. 

7.  Nondiffuse  scattering  is  most  likely  to  occur 
from  surfaces  that  are  smooth  and  free  of  gross 
contamination . 

8.  The  dispersion  of  nondiffuse  scattering  is 
greater  for  the  in-plane  pattern  than  for  the 
out-of-plane  pattern.  (That  is,  the  probability 
is  low  that  molecules  will  be  scattered  out  of 
the  plane  defined  by  the  incident  beam  and  the 
target  normal.) 

Notice  that  these  eight  characteristics  are  concerned 
with  gross  surface  properties  and  no  mention  is  made  of 
reflected  velocities  or  energy.  This  limitation  was  due  to 
the  unavailability  of  reliable  experimental  results  on  these 
quantities  at  that  time  (1967). 

The  reciprocity  condition,  in  general,  places  no 

restriction  upon  the  variations  of  reflected  characteristics 

with  changes  in  wall  or  gas  properties.  Thus,  no  definite 

statement  can  be  made  concerning  characteristics  2,  4,  5  and 

7.  However,  if  it  is  assumed  that  the  wall  properties  are 

independent  of  T  (say  a  and  a  are  constant)  then  a  limited 

w  x  y 

statement  may  be  made  concerning  characteristic  2.  The 

effect  of  T  is  manifest  in  S’.  As  T  increases,  S'  de- 
w  w 

creases  and  from  an  inspection  of  Figures  34,  35  and  36, 
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Appendix  A,  it  can  be  seen  that  characteristic  2  is  satisfied 
if  ax  1  cy 

Nothing  can  be  said  about  characteristics  4,  5  and  7 
based  upon  reciprocity  since  it  has  been  experimentally  deter¬ 
mined  that  the  accommodations  of  energy  between  the  gas  and 
surface  is  a  strong  function  of  m  (e.g.,  see  [33],  Figure  2). 

Based  upon  the  results  presented  in  Figures  34,  35 
and  36,  Appendix  A,  it  can  be  seen  that  the  present  theory 

satisfies  characteristic  1  if  a  £  a  .  Also  considering 

x  y 

Figure  37,  Appendix  A,  it  is  seen  that  characteristic  3  ap¬ 
pears  to  be  fulfilled  for  all  a  and  a  .  Likewise,  charac- 

x  y 

teristic  6  is  valid  if  a  <  a  from  an  investigation  of 

x  y 

Figure  34,  Appendix  A.  And,  finally,  characteristic  8  is 

confirmed  for  a  <  a  from  a  consideration  of  Figure  28, 
x  y 

Appendix  A. 

Thus,  the  present  gas-surface  interaction  model  is 
in  agreement  with  general  trends  as  presented  in  character¬ 
istics  1,  2,  3,  6  and  8  if  a  <  a  and  is  not  in  disagree- 

x  y 

ment  with  characteristics  4,  5  and  7. 

4.2  DETAILED  COMPARISON  WITH  EXPERIMENTAL  RESULTS 

4.2.1  Comparison  with  Data  of  Moran  [34] 

In  order  to  adequately  check  the  theoretical  results, 
it  is  desirable  to  have  experimental  measurements  on  several 
reflected  molecular  properties  from  a  given  surface  for 
various  incident  angles  and  velocities.  The  most  complete 
experimental  investigation  of  the  interaction  between  a 
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monatomic  gas  and  an  isotropic  surface  is  that  done  by 

Moran  [34]  for  argon  on  platinum.  He  measured  the  reflected 

time  of  flight  signal  as  a  function  of  reflected  angle  for 

two  incident  energies  at  several  incident  angles  from  which 

several  reflected  properties  are  determined.  His  results 

are  presented  in  Figures  38  through  44,  Appendix  A,  for 

reflected  relative  flux,  speed  ratio,  and  mean  speed,  from 

a  heated  platinum  surface. 

The  reflected  speed  ratio,  Sr,  is  a  measure  of  the 

thermal  spread  in  the  reemitted  molecular  velocity  distri- 

—  2 

bution.  Moran  defines  Sr  in  terms  of  v  and  v  as  follows: 

Sr  =  (^/v2-!)172  (209) 

Notice  that  any  consistent  error  in  the  scale  of  the 

velocities  tends  to  cancel  out  in  this  definition  of  S  . 

r 

As  has  been  pointed  out  previously  in  this  work,  the 
solution  for  n  =  0  in  rectangular  coordinators  (Equations 
13,  38  and  47)  seems  to  be  adequate  to  describe  most  experi¬ 
mental  results.  Values  of  a  and  a  were  chosen  so  as  to 

x  y 

best  fit  the  reflected  flux  distributions  in  a  least  squares 
sense.  Moran,  like  most  other  experimentalists,  presents 
flux  in  only  relative  terms,  i.e.,  arbitrary  units.  This  is 
due  to  the  difficult  problem  of  the  absolute  calibration 
of  a  molecular  flux  instrument.  The  relative  molecular  flux, 

'X. 

I',  is  calculated  from  I  by  introducing  an  aribtrary  multi¬ 
plication  factor  so  that  the  mean  value  of  I*  from  experiment 
and  from  theory  is  identical. 
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The  values  of  a  and  a  determined  for  the  best  fit 
x  y 

was  0.147  and  0.710,  respectively.  Using  these  values,  I', 

— 

v  and  v  can  be  determined  by  the  numerical  integrations  as 
discussed  in  Section  3.2,  and  Sr  obtained  from  Equation  209. 
These  theoretical  results  are  shown  as  the  curves  in  Figures 
38  through  44,  Appendix  A. 

4.2.2  Comparison  with  Data  of  Romney  [35] 

The  data  of  Moran  was  limited  to  low  incident 
velocity  (S'  =  0.742  and  1.313).  However,  Romney,  using  a 
seeded  beam  technique,  wherein  the  argon  was  accelerated 
using  a  light  carrier  gas,  obtains  reflected  density  distri¬ 
butions  of  high  velocity  argon  (S'  =  1.14  to  9.38)  reflected 
from  a  crystalline  silver  surface.  The  silver  surface  he 
used  was  vapor  deposited  in  situ  during  the  scattering 
measurements  so  as  to  ensure  an  uncontaminated  surface. 

Using  the  same  method  as  described  in  the  previous 

section,  a  and  a  were  determined  to  be  0.77  and  0.85, 
x  y 

respectively.  The  theoretical  results  for  the  relative 

reflected  density  distribution  for  an  incident  argon  beam 

with  energy  of  2.56  ev  is  compared  to  the  experimental 

result  in  Figure  45-a,  Appendix  A. 

Romney  also  determines  the  variation  of  the  position 

of  maximum  reflected  density  with  incident  energy  for  various 

angles  of  incident.  This  maximum  density,  for  the  conditions 

of  Romney's  test,  is  indistinguishable  from  the  position  of 

maximum  reflected  flux,  9,  therefore  his  results  will  be 

m 
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assumed  to  be  equal  to  0^.  Figure  46,  Appendix  A,  presents 

these  experimental  values  of  0m  for  incident  angles  of  35°, 

50°  and  65°  and  the  theoretical  results  obtained  using  the 

previously  determined  values  of  a  and  a  . 

x  y 

One  of  the  major  disadvantages  of  the  seeded  beam 
techniques,  as  used  by  Romney,  is  that  the  molecular  flux 
of  the  species  of  interest  is  drastically  reduced  for  the 
higher  energies,  thereby  making  the  measurement  of  reflected 
velocities  very  difficult.  At  the  present  there  are  no 
known  data  available  on  the  velocity  of  reflected  molecules 
under  high  energy  conditions. 

4.2.3  Comparison  with  Data  of  Bishara  [26,  36] 

Bishara  [29,  36]  has  obtained  data  on  the  flux  and 
mean  velocity  of  argon  scattered  from  a  crystalline  silver 
surface.  These  results  for  an  incident  energy  of  0.64  ev 
and  a  surface  temperature  of  550°K  are  shown  in  Figures  46, 
47  and  48,  Appendix  A,  for  incident  angles  of  -40°,  -50° 
and  -60°,  respectively.  He  made  measurements  not  only  in 
the  principal  plane  (qp  =  0)  but  at  two  transverse  planes 
(0=0  and  9  =  0  ) .  While  both  Romney  [35]  and  Bishara 
used  silver  surfaces,  Bishara  did  not  deposit  silver  in  situ 
as  did  Romney  but  polished  and  etched  the  surface  in  the 
atmosphere  prior  to  testing.  Data,  in  fact,  indicated  that 
the  scattering  pattern  varied  during  the  period  of  testing. 

Using  the  technique  previously  described  whereby  ax 
and  oty  were  determined  to  best  fit  the  experimental  data, 
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values  of  0.92  and  1.0  were  obtained,  respectively,  from 
Bishara’ s  results.  However,  the  theoretical  results  did 
not  adequately  explain  the  experimentally  observed  reflected 
flux  distributions.  Near  the  specular  direction  the  data 
were  more  peaked  while  further  away  the  data  seemed  to  be 
more  diffuse.  This  characteristic  is  similar  to  that  ob¬ 
served  by  Saltsburg  and  Smith  [37]  for  a  clean  silver  surface 
purposely  contaminated  by  exposure  to  the  atmosphere.  This 
explanation  is  consistent  with  the  higher  a's  obtained  from 
the  Bishara  [29]  data  than  the  Romney  [35]  data. 

For  any  real  surface,  especially  a  contaminated  one, 
there  are  variations  in  properties  over  the  surface,  i.e., 
the  surface  is  not  homogeneous.  Recalling  that  the  reci¬ 
procity  relation  is  linear,  the  scattering  kernel  could  be 
made  up  of  several  kernels  each  of  which  would  correspond 

to  a  micro-surface  and  have  an  a  and  a  associate  with  it. 

x  y 

In  regard  to  contaminated  surfaces,  recent  experi¬ 
mental  results  from  Dawbarn  et  al  [38]  for  nitrogen  on  con¬ 
taminated  aluminum  indicate  that  the  scattering  is  toward 
the  backscatter  direction  (i.e.,  Equation  5  is  applicable  to 
the  scattering  kernel). 

In  an  attempt  to  better  describe  theoretically  the 
results  of  Bishara  [29,  36]  it  was  assumed  that  the  scatter¬ 
ing  kernel  was  a  linear  combination  of  two  kernels  each  of 
which  could  be  reflected  according  to  Equation  5.  The  best 
fit  to  the  reflected  flux  distribution,  both  in  the  principal 
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plane  and  in  the  transverse  planes,  was  obtained  for  the 
following  kernel. 

P(v\v)  =  0.755  P[  +  0.245  Pg  (210) 

where: 

Pi  is  a  backscattered  kernel  with  a  =  0.99  and 

X  A 

Oy.  =  0.99 

P„  is  a  scattering  kernel  with  a  =  0.53  and  a  =  1.0 
a  x  y 

Using  Equation  210  the  mean  reflected  velocity  can  also  be 
obtained  by  numerical  integration,  these  results  for  v  and 
I'  are  compared  with  the  experimental  data  in  Figures  46 
through  48,  Appendix  A. 
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V.  SUMMARY  AND  CONCLUSIONS 

In  the  paper  by  Kuscer  et  al  [7]  the  statement  is 
made,  ’’It  is  surprisingly  difficult  to  construct  non-trivial 
mathematical  models  for  P(  v"' ,"v)  satisfying  all  three  require 
ments,"  (Equations  1,  2  and  4).  In  the  present  work,  by 
assuming  a  product  solution  for  P(T’ ,"v)  in  each  of  three 
coordinate  systems,  several  solutions  are  found  that  satisfy 
Equations  1,  2  and  4.  It  is  shown  that  one  particular  solu¬ 
tion,  as  given  by  Equations  38  and  47,  agrees  with  experi¬ 
mental  results  in  regard  to  reflected  flux  distributions, 
mean  velocities,  and  speed  ratio.  This  is  the  first  time 
that  any  theoretical  result  has  been  compared  to  all  these 
detailed  reflected  quantities.  The  theoretical  results  for 
the  position  of  the  maximum  reflected  flux  are  also  shown 
to  be  in  agreement  with  experimental  data,  and,  for  typical 
surfaces  with  high  incident  velocity,  the  reflected  lobe  is 
supra specular. 

A  more  general  scattering  kernel  can  be  obtained  by 
a  linear  combination  of  several,  or  in  the  limit  an  infinite 
number  of  individual  scattering  kernels.  In  the  limit  the 
infinite  sum  becomes  an  integral  and  the  following  can  be 
obtained  in  rectangular  coordinates: 


p<v’,v)"a  l  L  F(cvW  W'W  vv-w 

000  3  3  3  3 


W'W  dax  day  daz 


(211) 
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where  F  is  an  arbitrary  function  subject  to: 

111 

/  /  /  F(a,  a,  a)  da  da  da_  =  1  (212) 

000  x.  y  *  x  y  z 

and  such  that  P(v*,v)  is  non-negative.  Notice  that  this 
scattering  kernel  is  not,  in  general,  of  a  product  form 
unless  F  is  assumed  to  be  a  product.  The  scattering  kernel 
given  by  Equation  211  is  so  general  that  is  is  not  very 
helpful  for  the  description  of  a  practical  gas-surface 
interaction.  A  result  analogous  to  Equation  211  can  be 
obtained  in  cylindrical  and  spherical  coordinates,  and  the 
scattering  kernel  can  be  further  generalized  by  a  linear 
combination  of  the  results  in  various  coordinate  systems. 

Thus,  the  results  obtained  in  this  paper  when  com¬ 
bined,  using  the  linearity  and  reflectivity  (Equation  5) 
property,  permit  an  almost  unlimited  generalization.  How¬ 
ever,  reciprocity  places  restrictions  upon  the  variation  of 
the  scattering  kernel  with  incident  gas  conditions,  e.g., 

F  in  Equation  211  cannot  be  a  function  of  incident  gas 
energy  or  direction. 

The  primary  usefulness  of  the  results  obtained  herein 
is  not  the  "curve  fittings"  of  existing  data  per  se,  but  the 
determination  of  gas-surface  effects  under  conditions  where 
experimental  data  are  not  available.  In  other  words,  reci¬ 
procity  helps  us  extrapolate.  Thus,  surface  parameters  such 

as  a  and  a  obtained  under  moderate  incident  energy  condi- 
x  y 

tions,  can  be  used  to  obtain  gas-surface  interaction 
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parameters  under  high  energy  conditions. 

In  regard  to  future  experimental  investigations, 
it  is  desirable  that  investigators  publish  more  detailed 
results  on  the  reflected  molecular  distribution.  At  present, 
reflected  mean  molecular  quantities  are  determined  by  means 
of  a  numerical  integration,  with  appropriate  weighting 
factors,  of  the  time  of  flight  signal,  which  is  related  to 
P(T'  ,"v)  .  Much  of  the  details  of  the  gas-surface  interaction 
are  thus  lost  in  this  integration.  Experimental  results, 
e.g.,  [34,  35,  36],  in  the  past  have  usually  been  presented 
as  mean  reflected  quantities  mainly  because  of  the  unavail¬ 
ability  of  theories  which  could  be  compared  to  the  more 
detailed  results.  Therefore,  in  order  to  better  confirm  the 
present  theory  it  would  be  desirable  to  compare  direct  with 
the  experimental  scattering  kernel  determined  from  the  time 
of  flight  signal  of  the  reflected  molecules. 

As  more  experimental  data  are  obtained  under  more 
extreme  conditions  of  gas  and  surface,  it  would  be  necessary 
to  substantiate  the  applicability  of  Equations  38  and  47. 

It  may  be  that  under  some  conditions  that  other  scattering 
kernels,  developed  in  the  present  work,  are  necessary  to 
completely  describe  the  gas-surface  interactions. 
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Figure  1.  Coordinate  Systems. 
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Figure  2.  Plot  of  \/ 2RTW  Px>0. 
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Figure  2.  (Continued). 
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Figure  3.  (Continued). 
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Figure  4.  (Continued). 
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Figure  4.  (Continued). 
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Figure  5.  (Continued). 
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Figure  6.  (Continued) 


Figure  7.  (Continued). 
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Figure  7.  (Continued) 
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Figure  8.  (Continued). 
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Figure  8.  (Continued). 
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Figure  9.  (Continued). 


Figure  9.  (Continued). 
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Figure  9.  (Continued). 
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Figure  10.  Plot  of  2RTW  Pr  2. 
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Figure  10.  (Continued) 


AEDC -TR-72-158 


Figure  10.  (Continued) 
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Figure  11.  Plot  of  (2RTW)3'2  P, 
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Figure  12.  (Continued). 
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Figure  13.  (Continued). 
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Figure  18.  Sy  for  n  =  1  in  Spherical  Coordinates. 
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Figure  21.  %/§$  and  Vs^for  n  =  0  in  Rectangular  and  Cylindrical  Coordinates. 
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Figure  23.  n/s$  and  Vs?  for  n  =  2  in  Rectangular  and  Cylindrical  Coordinates. 


AEDC-TR-72-158 


Figure  25.  n/s*  for  n  =  1  in  Spherical  Coordinates. 
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Figure  27.  Molecular  Beam  and  Detector  Geometry. 


AEOC-TR-72-158 


AEDC-TR-72-158 


a.  <p  =  0,  0j  =  -20° 

Figure  28.  Reflected  Molecular  Intensity  (a*  =  0.5,  ay  =  0.8). 
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Figure  28.  (Continued). 
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c.  0  =  0,  0  =  -50°. 
Figure  28.  (Continued). 
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d.  0  =  8m,  6\  =  -50°. 
Figure  28.  (Continued). 
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e.  0  =  0,  0j  =  -80°. 
Figure  28.  (Continued). 
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f.  6-em,6i-S(f. 

Figure  28.  (Continued). 
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a.  <t>  =  0,  e,  =  -20°. 

Figure  29.  Reflected  Molecular  Intensity  (ax  =  0.8,  av  =  0.5). 
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b.  0  =  Om,  6,  =  -20°. 
Figure  29.  (Continued). 
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C.  0  =  0,  9  =  -50°. 
Figure  29.  (Continued). 
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e.  <p  =  0,di  =  -80°. 
Figure  29.  (Continued). 
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f-  Q  =  em.  dt  =  -so0. 
Figure  29.  (Continued). 
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a.  0  =  0,  0j  =  -20°. 

Figure  30.  Reflected  Molecular  Intensity  (ax  =  Oy  =  A-2). 
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b.  9  =  Qm,  0\  =  -20°. 
Figure  30.  (Continued). 
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c.  0  =  0 ,6  =  ■  50°. 
Figure  30.  (Continued). 
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(Continued). 
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Figure  30.  (Continued). 
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a.  <j>  =  0,  0j  =  -20°. 

Figure  31.  Mean  Reflected  Dimensionless  Speed  (a*  =  0.5,  ay  =  0.8). 
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b.  0=em,0 i  =  -20°. 
Figure  31.  (Continued). 
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c.  <p  =  0,  o,  =  -50°. 
Figure  31.  (Continued). 
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e.  0  =  0,  6;  =  -80°. 
Figure  31.  (Continued). 
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Figure  31.  (Continued). 
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Figure  32.  (Continued). 
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Figure  32.  (Continued). 


189 


AEDC-TR- 72-158 


Curve 

S’ 

0 

m 

X 

0.0 

41.78° 

2 

0.5 

48.60° 

3 

1 .0 

46.22° 

4 

1.5 

42.85° 

5 

2.0 

40.85° 

6 

2.5 

39.73° 

7 

3.0 

39.03° 

0.0  I _ I _ I _ 1 _ I _ I _ I 

-90°  -60°  -30°  0°  30°  60°  90° 

<F 


d.  9  =  0m ,  0,  =  -50°. 
Figure  32.  (Continued). 
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Figure  32.  (Continued). 


191 


AEDC-TR-72-1 58 


3.0 


2.0 


S 


1.0 


Curve 

S' 

e 

m 

1 

0.0 

.u 

oo 

o 

2 

0.5 

54.76° 

3 

1.0 

60.71° 

4 

1.5 

64.48° 

5 

2.0 

67.07° 

6 

2.5 

68.94° 

7 

3.0 

70.30° 

'  I _ I _ I _ I _ I _ I _ l 

-90°  -60°  -30°  0°  30°  60°  90° 

9 


f.  6  =  6m,  6-,  =  -  80°. 
Figure  32.  {Continued). 
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0 


a.  <p  =  0,  6,  =  -20°. 

Figure  33.  Reflected  Value  of  AS  (a*  =  Oy  =  A  2). 
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b.  e  =  em,ei  =  . 20°. 

Figure  33.  (Continued). 
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c.  4>  =  0,6i  =  -  50°. 

Figure  33.  (Continued). 
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d.  6  =  0m,0i  =  -50°. 
Figure  33.  (Continued). 
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e.  0  =  0,  0;  =  -80°. 
Figure  33.  (Continued). 
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f.  e  =  em,e j  =  -80°. 

Figure  33.  (Continued). 
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a.  Relative  Flux. 

Figure  38.  Comparison  of  Theory  and  Data  for  Argon  Scattered  from 
Platinum,  0-,  =  -40°,  v'  =  5.4  x  IQ4  cm/sec. 
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b.  Speed  Ratio. 


c.  Mean  Speed. 
Figure  38.  (Continued). 
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a.  Relative  Flux. 

Figure  39.  Comparison  of  Theory  and  Data  for  Argon  Scattered  from 
Platinum,  Q\  =  -55°,  v'  =  5.4  x  104  cm/sec. 
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d.  Speed  Ratio. 


c.  Mean  Speed. 
Figure  39.  (Continued). 
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a.  Relative  Flux. 

Figure  40.  Comparison  of  Theory  and  Data  for  Argon  Scattered  from 
Platium,  0j  =  -65°,  v'  =  5.4  x  104  cm/sec. 
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b.  Speed  Ratio. 


e 


c.  Mean  Speed. 
Figure  40  (Continued). 
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a.  Relative  Flux. 

Figure  41.  Comparison  of  Theory  and  Data  for  Argon  Scattered  from 
Platinum,  Q\  =  -20° ,  v'  =  9.55  x  104  cm/sec. 
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b.  Speed  Ratio. 


c.  Mean  Speed. 
Figure  41.  (Continued). 
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a.  Relative  Flux. 

Figure  42.  Comparison  of  Theory  and  Data  for  Argon  Scattered  from 
Platinum,  0j  =  -40°,  v'  =  9.55  x  104  cm/sec. 
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b.  Speed  Ratio. 


c.  Mean  Speed. 
Figure  42.  (Continued). 
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a.  Relative  Flux. 

Figure  43.  Comparison  of  Theory  and  Data  for  Argon  Scattered  from 
Platinum,  0S  =  -65°,  v'  =  9.55  x  104  cm/sec. 


213 


AEDC-TR-72-1  58 


O  Data,  Koran  [34] 


b.  Speed  Ratio 


c.  Mean  Speed. 
Figure  43.  (Continued). 
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a.  Relative  Flux. 

Figure  44.  Comparison  of  Theory  and  Data  for  Argon  Scattered  from 
Platinum,  =  -80°,  v'  =  9.55  x  104  cm/sec. 
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b.  Speed  Ratio. 


c.  Mean  Speed. 
Figure  44.  (Continued). 
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Data,  Romney  [35] 
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a.  Relative  Density. 

Figure  45.  Comparison  of  Theory  with  Data  of  Romney  [35] . 


b.  Position  of  Maximum  intensity. 
Figure  45.  (Continued). 
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c.  Relative  Intensity,  0  =  0°. 
Figure  46.  (Continued). 
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Data,  Bishara  [29,35] 
Tw  -  550°K 

Theory,  Based  Upon 
Equation  (210) 


O  Data,  Bisbara  [29,35] 
Tw  -  550°K 

- Theory,  Based  Upon 


c.  Relative  Intensity,  0-0°. 
Figure  47.  (Continued). 
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d.  Mean  Speed,  0  =  0°. 


c.  Relative  Intensity,  0  =  0°. 
Figure  48.  (Continued). 
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APPENDIX  B 

EVALUATION  OF  CERTAIN  INTEGRALS  OF  EXPONENTIALS 
AND  BESSEL  FUNCTIONS 

At  several  points  in  the  development  of  the  analysis 
it  was  necessary  to  evaluate  infinite  integrals  of  the 
following  type: 

00  n  O 

/  xn  exp(-ax^)  iy(bx)  dx  (B-l) 

0 

1/2 

Letting  x  =  t  expression  B-l  becomes: 

1/2  /  tn/2-1/2  exp(-at)  Iv(bt1/2)  dt  (B-2) 

0 

This  integral  can  be  evaluated  by  the  use  of  Equation 
6.643-2  of  Gradshteyn  and  Ryzhik  [39]  and  the  result  is: 

r(n/2  +  v/2  +  1/2)  [an/2  b  r(v+l)]_1  exp(b2/8a)M_n/2j  v/2(b2/4a) 

(B-3) 

where  M  .  (t)  is  Whittaker's  Function  (see  [30],  page  505). 

The  result  in  expression  B-3  can  be  written  in  terms 
of  Rummer's  Function  M(c,d,t)  using  Equation  13.1.32  of  [30] 
and  the  final  result  for  the  integral  is: 

r(n/2  +  v/2  +  1/2)  [an/2  b  r(v+l)]_1(b2/4a)  (1't'V)/2  • 
M(l/2  +  v/2  +  n/2,  1  +  v,  b2/4a)  (B-4) 
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where  from  [30]  Rummer's  Function  is  defined  as: 


H(e.d,t)  =  V  r .(g5?_r.M_g 
r(c)  r(d+m)  m.' 

m=0 


Notice  that  for  c  =  d  Equation  B-5  becomes; 


M(d,d,t)  -  exp(t) 

Using  these  results  the  following  equations  are 
obtained  for  some  integrals  of  interest; 


00  2  n 

/  xe”ax  IQ(bx)dx  -  l/(2a)  exp(bZ/4a) 


00  2 

I  x2e"ax  I0(bx)  dx  =  7r1/2/(4a3/2>  M(3/2, 1,  b2/4a) 

0 

/  xe“ax2  I^bx)  dx  -  (b7T1/2/8a3/2)  M(3/2, 2, b2/4a) 

0 

00  n  2  _  _ 
f  x  e”ax  I^(bx)  dx  =  b/(4aZ)  exp(bZ/4a) 

0 

/  x3e"ax2  ri(bx)  dx  =  (37r1/2b/16a5/2)  M(5/2, 2,  b2/4a) 
0 


(B-5) 

(B-6) 

(B-7) 

(B-8) 

(B-9) 

(B-10) 

(B-ll) 
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